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ABSTRACT 

We have investigated the observational characteristics of a class of broad emission line 
region (BLR) geometries that connect the outer accretion disc with the inner edge 
of the dusty toroidal obscuring region (TOR). We suggest that the BLR consists of 
photoionised gas of densities which allow for efficient cooling by UV/optical emission 
lines and of incident continuum fluxes which discourage the formation of grains, and 
that such gas occupies the range of distance and scale height between the continuum- 
emitting accretion disc and the dusty TOR. As a first approximation, we assume a 
population of clouds illuminated by ionising photons from the central source, with the 
scale height of the illuminated clouds growing with increasing radial distance, forming 

■ an effective surface of a " bowl" . Observer lines of sight which peer into the bowl lead 

to a Type 1 Active Galactic Nuclei (AGN) spectrum. We assume the gas dynamics are 

q , dominated by gravity, and we include in this model the effects of transverse Doppler 

■ shift, gravitational redshift and scale- height dependent macro-turbulence. 

^ \ Our simple model reproduces many of the commonly observed phenomena as- 

sociated with the central regions of AGN, including : (i) the shorter than expected 
continuum-dust delays (geometry), (ii) the absence of response in the core of the 
optical recombination lines on short timescales (geometry /photoionisation), (hi) an 
enhanced red-wing response on short timescales (GR and TDS), (iv) the observed 
| differences between the delays for high- and low-ionisation lines (photoionisation), (v) 

£C) , identifying one of the possible primary contributors to the observed line widths for 

' near face-on systems even for purely transverse motion (GR and TDS), (vi) a mecha- 

\ nism responsible for producing Lorentzian profiles (especially in the Balmer and Mg II 

emission-lines) in low inclination systems (turbulence), (vii) the absence of significant 
£^ ' continuum-emission-line delays between the line wings and line core (turbulence; such 

£v^| \ time-delays are weak for virialised motion, and turbulence serves to reduce any dif- 

ferences which may be present), (viii) associating the boundary between population 
A and population B sources (Sulentic et al. 2000) as the cross-over between incli- 
nation dependent (population A) and inclination independent (population B) line 

■ profiles (GR+TDS), (ix) provides a partial explanation of the differences between the 
' emission-line profiles, here explained in terms of their line formation radius (photoion- 

isation and/or turbulence), and (x) the unexpectedly high (but necessary) covering 
fractions (geometry). 

A key motivation of this work was to reveal the physical underpinnings of the 
reported measurements of supermassive black hole (SMBH) masses and their uncer- 
tainties. We have driven our model with simulated continuum light-curves in order to 
determine the virial scale factor /, from measurements of the simulated continuum- 
emission-line delay, and the width (fwhm, 07) and shape (fwhm/cri) of the rms and 
mean line profiles for the energetically more important broad UV and optical recom- 
bination lines used in SMBH mass determinations. We thus attempt to illuminate the 
physical dependencies of the empirically determined value of /. We find that SMBH 
masses derived from measurements of the fwhm of the mean and rms profiles show 
the closest correspondence between the emission lines in a single object, even though 
the emission line fwhm is a more biased mass indicator with respect to inclination. 
The predicted large discrepancies in the SMBH mass estimates between emission lines 
at low inclination, as derived using 07, we suggest may be used as a means of iden- 
tifying near face-on systems. Our general results do not depend on specific choices 
in the simplifying assumptions, but are in fact generic properties of BLR geometries 
with axial symmetry that span a substantial range in radially-increasing scale height 
supported by turbulence, which then merge into the inner dusty TOR. 

Key words: methods : numerical - line : profiles - galaxies : active - quasars : 
emission lines 
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1 INTRODUCTION 

Intensive multi- wavelength monitoring campaigns of a hand- 
ful of individual AGN have radically altered our view of 
the broad emission- line region (hereafter, BLR). Correlated 
multi-wavelength continuum and emission-line variations 
(reverberation mapping, hereafter RM) reveal that the BLR 
is small in size (i^BLR oc L^'v), and shows strong gradients in 
density and/or ionisation parameter (e.g., Clavel et al. 1991; 
Peterson et al. 1991; Krolik et al. 1991). The gas dynamics 
are largely dominated by the central super-massive black 
hole, a realisation which when coupled with estimates of the 
BLR size and gas velocity dispersion, have enabled the deter- 
mination of virial black hole masses in « 40 nearby AGN 
(Peterson 2010). The derived (geometry dependent) virial 
scale factors appear approximately constant amongst the 
emission lines in individual sources, and over many seasons 
for which the mean source luminosity can differ considerably 
(see Krolik et al. 1991; Peterson & Wandel 1999, 2000; Peter- 
son et al. 2004; Kollatschny 2003; Bentz et al. 2007; Denney 
et al. 2010; Peterson 2010, and references therein), consistent 
with a virialised velocity field within an ionisation stratified 
BLR. Further, the mass estimates compare favourably (to 
within a factor of a few) with those determined from stellar 
velocity dispersion estimates in nearby AGN (e.g. Ferrarese 
et al. 2001), and importantly allow scaling relations (e.g., 
BLR size vs. luminosity) to be derived enabling characteri- 
sation of black hole masses in AGN with only single epoch 
spectra (e.g., see Park et al. 2012). 

A key result of previous RM campaigns is that the 
recovered response functions for the optical recombination 
lines do not reach a maximum at zero delay as might be 
expected for a spherically symmetric distribution of line 
emitting gas (see e.g. Home, Welsh & Peterson 1991; Fer- 
land et al. 1992; Home, Korista & Goad 2003; Bentz et al. 
2010), which when taken at face- value suggests that the BLR 
gas lies well away from the line of sight to the observer, 
prompting a number of geometry dependent interpretations 
(for example, see Manucci, Salvati and Stanga 1992; Chi- 
ang and Murray 1996; Murray and Chiang 1997; Bottorff 
et al. 1997). Alternative explanations include optically-thick 
line emission (see e.g. Ferland, Shields, & Netzer 1979; Fer- 
land et al. 1992; Goad 1995, Ph.D. thesis; O'Brien, Goad 
& Gondhalekar 1994) or the presence of low-responsivity 
gas in the inner BLR (e.g. Sparke 1993; Goad, O'Brien & 
Gondhalekar 1993). To break the model degeneracy requires 
determination of ^(v,t) (Welsh & Home 1991; Goad et al. 
1995, Ph.D. thesis; Home et al. 2004). Initial attempts were 
limited to determining the "time- averaged response" as a 
function of several bins in line of sight velocity within the 
Balmer, He I A5876 and He n A4686 emission line profiles 
(see e.g., Ulrich & Home 1996; Kollatschny & Bischoff 2002; 
Kollatschny 2003; Kollatschny & Zetzl 2010; Denney et al. 
2009, 2010; Bentz et al. 2010a) with somewhat mixed re- 
sults. However, improvements in data quality and the con- 
tinued refinement of inversion techniques have now begun 
to reveal structure in the 2-d response functions of the opti- 
cal recombination lines in Arp 151 (Bentz et al. 2010b) and 
NGC 4051 (Denney et al. 2011). For Arp 151, V(v,t) shows 
that the red-wing responds first and, significantly, a deficit 
in response in the core of the lines at short time-delays. 
These data have been used to exclude models with radial 



outflow, as well as flattened discs and thick spherical shell 
geometries, and demonstrate consistency with a warped disc 
geometry. On the other hand, Denney et al. (2011) suggest 
that in NGC 4051, *Sf(y,r) for H/3 is consistent with a disc- 
like BLR geometry with the response weighted toward larger 
BLR radii, and possibly a scale-height dependent turbulent 
component ( N3.3p . 

1.1 The BLR outer radius 

At high incident continuum photon fluxes, dust grains can- 
not survive and atomic emission lines dominate the cooling 
of the gas. Conversely, at low incident fluxes heavy elements 
condense into grains which absorb much of the incident con- 
tinuum flux and become an important cooling agent of the 
gas, leading to a diminution of the line emission. Thus, as 
was first suggested by Netzer & Laor (1993), the outer edge 
of the BLR is likely largely determined by the distance at 
which grains can survive (see also Nenkova et al. 2008; Mor 
& Trakhtenbrot 2011; Mor & Netzer 2012), marking the in- 
ner boundary of the "toroidal obscuring region" (hereafter, 
TOR). Several studies have concentrated on determining the 
distance to the hot dust, from measurements of the delay 
between the UV/optical continuum and near-IR (—2 /mi) 
thermal emission (e.g., Minezaki et al. 2004; Suganuma et al. 
2004, 2006; Yoshii et al. 2004; Kishimoto et al. 2007; Koshida 
et al. 2009). These reveal that Tdust oc L 0-5 , and is uncorre- 
cted with black hole mass. However, surprisingly Tdust is a 
factor of —2-3 smaller than expected from the grain subli- 
mation radius (Barvainis 1987; Nenkova 2008; Suganuma et 
al. 2006; Kishimoto et al. 2007) based on the central source 
luminosity and the assumption that R su b ~ crdust- For ex- 
ample, in NGC 5548 crdust ~ 50 light-days (Suganuma et al. 
2004; 2006), to be compared with an expected sublimation 
radius of -150 light days (Nenkova et al. 2008; Eq. 1). When 
compared with the lag measurements for the broad C in] and 
Mg II emission-lines in this source 28 and 34-72 days (with 
significant uncertainties), respectively (Clavel et al. 1991), 
places the outer regions of the BLR in the vicinity of the 
inner edge of the dusty torus. 

Hu et al. (2008), Zhu et al. (2009), Shields et al. (2010), 
and Mor & Netzer (2012) suggest that the outer BLR and 
inner torus overlap in an "intermediate" line region (ILR). 
In addition, Landt et al. (2011;conf.proc), adopting an 
isotropic, virialised velocity field for the BLR and black holes 
masses inferred from reverberation campaigns, used their 
sample of type 1 AGN with IR spectra (Landt et al. 2008) 
to show that the outer edge of the BLR scales as L 0-5 , and 
corresponds to incident continuum photon fluxes that are 
consistent with that expected at or just interior to the hot 
dust radius. 

Recent observational and theoretical studies of the TOR 
(e.g., Honig et al. 2006; Tristam et al. 2007; Nenkova et al. 
2008; Ramos Almeida et al. 2009) find that it is most likely 
clumpy, consisting of dusty clouds. Given that there must 
exist a reservoir of gas feeding the central accretion disc, it 
is natural to assume that the BLR consists of gas with a 
range of density and incident continuum flux which max- 
imises atomic line emission over atomic continuum emission 
(the inner accretion disc) and grain emission (TOR). The 
BLR might then be expected to span between the inner ac- 
cretion disc and the TOR both in radius and in scale height. 
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1.2 Geometry-dependent time-delays 

If interpreted in terms of a characteristic "size" for the 
dusty torus, then such small delays pose severe problems 
both for grain survival, requiring both larger and more ro- 
bust grains (e.g., graphite; Mor & Trakhtenbrot 201 lj3 or 
significant shielding. Furthermore, the reduction in radius 
that this implies for the outer BLR poses a severe problem 
for (the already large) BLR gas covering fractions (see e.g. 
Kaspi & Netzer 1999 and Korista & Goad 2000) required by 
photoionisation models of the broad emission-lines in this 
source, the outer radius of which was chosen to be roughly 
coincident with the distance at which hot dust can survive. 

The simplest way of reconciling the dust-delays with 
the expected dust formation radius is to invoke a geome- 
try for the dust which departs from spherical symmetry as 
was recently proposed by Kawaguchi & Mori (2010, 2011; 
see also Liu & Zhang 2011). In their model anisotropic con- 
tinuum radiation (see e.g. Netzer 1987; but see Nemmen & 
Brotherton 2010) acts through grain sublimation to impose 
a strong polar angle dependence on the distance at which 
grains can form (smaller at large polar angles), carving the 
TOR into a bowl-shaped geometry of significant scale-height 
H/R x « 1, where R x is measured along the roughly coinci- 
dent mid-planes of the accretion disc and the dusty TOR. 
This places the dusty gas nearer to the line of sight for ob- 
servers with typical viewing angles of 10-40 degrees with 
respect to the polar axis (for observers of type 1 AGN) , and 
more distant from the central continuum source (R), than a 
comparable ring of dust of radius R x lying within the mid- 
plane (for the same time-delay producing dust-delays that 
are substantially smaller than R su b as computed from the 
observed continuum. Curiously, missing from their model is 
the BLR, which in the context of the unified model of AGN 
(Antonucci 1993; Urry & Padovani 1995), lies somewhere 
between the inner accretion disc responsible for the optical- 
UV- X-ray continuum (R < 100 R g ), where R g = GM B h/c 2 , 
and the dusty obscuring torus (R > 20,000 R g ). A BLR with 
a significant scale height would cut off much of the interior 
of their model dusty torus from view of the incident contin- 
uum, obviating the need to minimise contributions to the 
hot dust response from very short time delays. It would also 
relegate most of the hot dust emission to dusty gas having 
greater scale heights than and lying at distances somewhat 
further from the central light source than the outer BLR, 
with a covering fraction smaller than one might infer from 
the overall structure of the TOR. Landt et al. (2011b) and 
Mor & Trakhtenbrot (2011) deduced covering fractions (typ- 
ically < 10% and ~ 13%, respectively) for the gas emitting 
the thermal emission from the hot grains, smaller than ex- 
pected from type 1/type 2 statistics (Schmitt et al. 2001; 
Hao et al. 2005) and simple geometrical considerations. The 
hot grains might lie along the top rim of a bowl-like geom- 
etry of clouds constituting the BLR and the dusty TOR. 

Spatially-resolved near-IR imaging have largely con- 
firmed the above basic picture of the TOR (Nenkova et al. 



1 Even graphite grains with sublimation temperatures of £^1700K 
with typical sizes of 0.1 microns would have a sublimation dis- 
tance of « 110(L bo j/10 44 - 3 ) - 5 light-days (Nenkova et al. 2008). 

2 Such a ring would have an observed delay of R x /c for an ob- 
server with a line of sight near the polar axis. 



2008) - clumpy with significant scale height (Krolik & Begel- 
man 1988) and interferometric ring radii that are typically 
a factor of ^2 or so larger than crdust (Honig et al. 2010; 
Kishimoto et al. 2007,2009a,b,2011; Raban et al. 2009; Tris- 
tram et al. 2007, 2009). 

1.3 Profile shape as a discriminator of BLR 
geometry 

For nearby AGN a comparison between their RM black hole 
mass estimates and those derived from their stellar velocity 
dispersion, indicate significant differences between the virial 
scale factors of the two broad emission lines often used, be- 
ing smaller for Hf3 than for C IV (Decarli et al. 2008). Decarli 
et al. (2008) argue that the smaller virial scale factor for Civ 
indicates a more flattened distribution for the gas producing 
this line, while H/3 originates in a more spherical distribu- 
tion, or equivalently a region with larger scale height (e.g., a 
thick disc). Sulentic et al. (2000) suggested that the absence 
of a correlation between the line widths of C IV and H/3 
in individual sources, together with observed differences in 
their line profile shapes, fwhm/a, indicate that C IV is pro- 
duced in a more flattened configuration, while H/3 originates 
from a region with significant scale height. The absence of a 
correlation in their respective line widths, is then explained 
in terms of an additional scale height-dependent turbulent 
velocity component for H/3 (see also Collin et al. 2006). Sig- 
nificantly, recent dynamical modelling of the nearby Sy 1 
Mrk 50 by the Lick AGN Monitoring Project (LAMP), is 
suggestive of a thick disc geometry for the Hf3 line emitting 
region (Pancoast et al. 2012). 

Fine et al. (2008, 2010, 2011) provide further evidence 
against a pure-planar/disc-like geometry for the line emit- 
ting gas. By combining line velocity dispersion data from 
four large quasar spectroscopic surveys they found that the 
dispersion in the distribution of broad-line widths is smaller 
for C IV than Mg II, and unlike Mg n is only weakly cor- 
related with source luminosity. While differences in the line 
widths of these two lines are consistent with ionisation strat- 
ification within a virialised gas, the small dispersion in the 
distribution of the line widths of both lines argues against a 
pure-planar/disc- like geometry for the line emitting gas. 

In a separate study, Kollatschny and Zetzl (2011) mod- 
elled the broad H/3 and C IV emission-line profiles of AGN 
with RM data, by convolving a rotational velocity com- 
ponent representing the gravitational potential of the cen- 
tral super-massive black hole, with a Lorentzian component, 
which they associated with turbulence. While for a given 
line, the degree of turbulence remains similar from one ob- 
ject to the next, they found that their model fits require a 
larger turbulent velocity component for C IV than for H/3. 
Assuming that the magnitude of the turbulent component 
increases with increasing scale- height, this argues for a larger 
scale height for C IV, in stark contrast to the results above. 



1.4 Summary and a proposed geometrical model 
of the BLR 

There is substantive evidence that the BLR geometry is nei- 
ther very flat nor spherical, with predominantly virialised 
gas dynamics, and is largely absent from regions lying near 
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to our privileged line of sight down into the opening of the 
TOR in the unified model of AGN. Nenkova et al. (2008) 
suggest that what divides the TOR from the BLR is the abil- 
ity of grains to survive or not within clouds moving through 
the high radiation energy density environments found near 
the central engines of AGN. Following Nenkova et al. (2008), 
we propose a model in which the BLR bridges the gap be- 
tween the inner continuum-emitting accretion disc and the 
hot grain emitting inner TOR (Netzer & Laor 1993), both in 
spatial distance from the super-massive black hole and also 
in characteristic scale height that grows with increasing dis- 
tance. Recent results from dynamical mass estimates for the 
mass of the central black hole in the nearby Seyfert 1 galaxy 
Mrk 50 (Pancoast et al. 2012), appear to favour dynamical 
models in which the BLR has substantial scale-height (ie. a 
thick disc geometry). Together the three structures present 
an effective surface of a bowl-like geometry - roughly flat 
on the bottom and scale height increasing with distance R x . 
The inner accretion disc and BLR form the bottom and 
lower sides of the bowl, while the gas containing the hot 
dust is then found on the upper reaches of the bowl. A simi- 
lar configuration has also been suggested recently by Gaskell 
(2009). 

In this work we explore the observational characteris- 
tics of bowl-shaped BLR geometries comparing our detailed 
model calculations with observations taken from the liter- 
ature. The outline of this work is as follows. In section 2 
we present generic toy bowl-shaped BLR models, illustrat- 
ing the dependence of the steady-state 2-d and 1-d response 
functions and emission-line profiles on the key parameters 
describing the model. In 33 we outline a fiducial BLR ge- 
ometry the parameterisation of which has been chosen to 
to be broadly representative of the conditions within the 
well- studied Seyfert 1 galaxy NGC 5548. Our model in- 
cludes a treatment of the transverse Doppler shift, grav- 
itational redshift and scale-height dependent macroscopic 
turbulence. In we present the results from photoionisa- 
tion model calculations of the radial surface emissivities for 
four broad emission lines commonly used in measuring AGN 
black hole masses. We also describe the model continuum 
light-curves, chosen to match the continuum variability be- 
haviour of NGC 5548, and used to drive our fiducial BLR 
geometry to produce time- variable emission-line profiles and 
light-curves. The results from our simulations are presented 
in 33 In we compare the results of our simulations with 
observations presented in the literature. 



2 A BOWL-SHAPED BLR GEOMETRY 

In order to investigate the properties (gas distribution, kine- 
matics, 1-d and 2-d response, profile shape etc) of bowl- 
shaped geometries as a class, we begin by first constructing 
a very simple toy model. For the purposes of illustrating 
the concept, we initially set the inner radius at 1 lt-day and 
outer radius at 10 It- days and assume that the clouds lying 
near the surface of the bowl radiate isotropically, and do not 
cover one another (ie., have a direct view of the central con- 
tinuum source). We parameterise the shape of the bowl in 
terms of the BLR scale height H (in units of lt-days), such 
that 




inner accretion disc R 
(continuum cooling) 



Figure 1. An illustration of the Bowl-shaped geometry. The BLR 
clouds occupy a region where gas is strongly illuminated by the 
central continuum source and line cooling dominates over either 
continuum cooling and grain cooling, effectively bridging the re- 
gion between the outer accretion disc and the hot dust. The 
quoted hot dust covering fraction (f2/47r ~ 0.1) is taken from 
Landt et al. 2011b. 



H = f3R™, (1) 

where f5 and a are constants, i? x is the projected radial 
distance (in lt-days) along the plane perpendicular to the 
observers line of sight, ie. R x = R(sin 0), where is the polar 
angle (see figure [1]). In the limiting case of zero scale- height, 
the normalisation constant /3 is zero by definition, and the 
material lies along the mid-plane of the disc. For fixed (non- 
zero) /?, the shape of the bowl is determined by the power- 
law index a. The bowl is convex for a < 1, concave for 
a > 1, and cone-shaped for a = 1 (figure 2a). Note that (f) 
is a function of radial distance R unless a — 1.0, for which 
(f) = constant. When viewed along the axis of symmetry 
(down into the bowl), the spread in time-delays is minimised 
for fixed /?, when a = 2, i.e. parabolic. For a > 2, gas lying 
near the surface of the bowl at larger radii will respond 
before gas at small radii (figure 2b). Since the time-delay, 
r, at radial distance R is r = R — H, this condition will be 
met provided SH > 5R (or equivalently dr/dR is negative). 
Note that in our model, we only observe the inner surface 
of the bowl, the outer surface being largely obscured from 
sight (since such sight lines would have to pass through the 
dusty TOR, see e.g. figure 1). 

The toy models discussed above, and illustrated in fig- 
ure 2 serve as a tool for exploring how specific combinations 
of a, P determine the bowl shape, and the expected spread 
in time-delay for a fixed R ou t- Since the formulism described 
leads to an infinite number of bowl shaped geometries, we 
restrict the parameter space further, by fixing both the BLR 
inner and outer radius (Ri n , Rout), and the expected delay 
at the outer radius, r(R = R ou t) (for i = 0). Then, for a 
given choice of a, f3 is adjusted to match the specified de- 
lay at the outer boundary. In so doing, we can examine the 
properties of families of bowl-shaped geometries, with fixed 
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Figure 2. Bowl-shaped geometries with inner radius Ri n = 1 lt- 
day, outer radius Rout = 10 lt-days, inclination i = 0, and scale 
height H (lt-days) and indicating (i) the shape (upper panels), 
and (ii) the range in expected delays (lower panels), for fixed ft 
(figures 2a, 2b), and fixed a (figures 2c, 2d). 

radial extent, and with shapes differentiated by their chosen 
value of a. 

2.1 The BLR velocity field 

For the velocity field, we assume that the gas motion is domi- 
nated by the gravitational potential well in which it sits, and 
is largely circularised, as might be expected in the presence 
of significant dissipative forces, such that: 



Vke P - A (jR 2 + 2i2 2a)3/2 ' ^ 

where Vke P is the local Keplerian velocity, K = GMbh, and 
Mbh is the mass of the central black hole. For the limiting 
case of a geometrically thin disc (ie. ft = 0), v\ ev — Kj^, 
and the velocity field reduces to that expected for simple 
planar Keplerian orbits. Note that in this formulism, there 
is no radial component to the velocity field. Significant ra- 
dial motion, either in the form of turbulence, or bulk radial 
motion may be introduced by the addition of an azimuthal 
perturbation to the velocity field (see £13.3[) . 

While we do not exclude the possibility of infall, out- 
flow (either centrifugal or radiation pressure dominated), or 
line-driven wind contributions to the BLR kinematics within 
individual AGN, we consider here the simplest possible 



and likely underlying scenario, of a largely gravitationally- 
dominated circularised flow. In particular we do not here 
consider the effects of a substantial wind contribution to the 
broad emission-line profiles of especially the high ionisation 
resonance lines, such as C IV (see e.g., the disc- wind models 
of Chiang & Murray 2006; Murray & Chiang 2007), which 
may be especially important in these lines in high Edding- 
ton rate AGN. We also note that while outflows, as mani- 
fested in blue-shifted UV and X-ray absorption lines (e.g., 
Elvis 2000) are common in AGN, their origin (s) is presently 
uncertain, although the inner TOR remains a plausible can- 
didate (Pier & Voit 1995; Mullaney et al. 2008). Indeed the 
velocities reported for these outflows in Seyfert AGN (up 
to ~1000 km s _1 , e.g., Crenshaw et al. 1999) roughly cor- 
respond with the expected gravitationally dominated veloc- 
ities at the purported distance of the TOR. Moreover, in 
the sample of AGN with reverberation mapping data (e.g., 
Kaspi et al. 2000; Richards et al. 2011), there is no evidence 
for strong blue asymmetries in their C IV emission-line pro- 
files, normally taken to be a strong indicator for outflows 
within the BLR. 

Rather, we introduce macro-turbulent cloud motion 
(§3.3) as a simple, and often suggested, mechanism for pro- 
viding a substantial radially-dependent scale height for the 
BLR clouds, allowing the BLR gas to intercept a signifi- 
cant fraction of the ionising continuum radiation necessary 
to produce the observed line strengths, something which re- 
mains problematic for thin-disc models of the BLR. In jus- 
tifying our assumption, we note that the TOR which is the 
likely reservoir of gas that ultimately feeds the accretion 
disc, itself has a very large scale- height, while dynamical 
models of the BLR based on reverberation data appear to 
favour flared-disc geometries with substantial opening an- 
gles above the disc mid-plane (e.g. Pancoast et al. 2012, see 
§3.3 for details). Importantly in 33.21 we do investigate the 
effects of both transverse Doppler shift and gravitational red- 
shift on the observed emission line profile and velocity delay 
map. A particular choice of Mbh, Rin, Rout, a, and r(R ou t) 
completely describes the BLR geometry and velocity field. 

We emphasise that we do not advocate a fixed inner and 
outer boundary for the BLR nor do we mean to imply that 
the BLR forms the smooth inner surface of a bowl-shaped 
geometry. These are merely the simplest parameterisations 
of something which is very likely far more complex0 Rather 
we aim to explore the observational consequences of assum- 
ing a bowl shaped BLR geometry in which the gas dynamics 
are dominated by the central super-massive black hole, and 
how these then impact on our interpretation of line profile 



3 The location of the effective (line emissivity weighted) inner 
and outer BLR radii will vary in response to the incident con- 
tinuum flux, being generally forced to larger radii in higher con- 
tinuum flux states. The inner radius is set by some combination 
of over-ionisation and emission line optical depth effects, along 
with line-width visibility in the deep potential well, presuming 
the availability of emitting gas. These adjustment time scales 
should be fairly rapid compared to the central continuum vari- 
ability time scales. The outer radius is set by the effects of grain 
heating/cooling and line destruction on emission line emissivity, 
and the grain vaporisation and condensation time scales which 
are likely comparable to or longer than, typical incident central 
continuum variability time scales. 
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shapes, correlated continuum and emission-line variability, 
velocity resolved response functions, black hole mass esti- 
mates and virial scale factors reported in the literature for 
both individual sources and among the AGN population as 
a whole. 



3 A FIDUCIAL BLR GEOMETRY 

To illustrate the general properties of bowl-shaped BLR ge- 
ometries we construct a fiducial bowl-shaped geometry for 
the BLR for comparison with observations. For expedience 
we choose a parameter set appropriate for the Seyfert 1 
galaxy NGC 5548. Our fiducial BLR geometry has a cen- 
tral black hole mass Mbh = 1.0 x 10 8 M©. We set the 
inner radius at 200 gravitational radii (^1.14 light-days for 
our assumed black hole mass) noting that the response time 
scales for He n A1640A and N V A1240A are ^2 days in 
NGC 5548 (Korista et al. 1995). We choose an outer BLR 
radius marking the upper rim of our bowl geometry of 100 
light-days, roughly speaking the graphite grain sublimation 
radius (Nenkova et al. 2008) for our chosen continuum nor- 
malisation (see £|4.ip , and a maximum time-delay at the 
outer radius (for i = 0) of r = (R — H)/c — 50 days, similar 
to the dust reverberation time-delay measured for this ob- 
ject. With these parameters, the source covering fraction as 
determined from the polar angle to the bowl rim (60 degrees) 
is 50% for our fiducial BLR geometry. 

We populate the bowl surface with discrete line emit- 
ting entities (hereafter, clouds) of fixed column density so 
that each cloud has an un-obscured view of the continuum 
source and radiates energy in a manner approximating the 
phases of the moon (see e.g., O'Brien et al. 1995, Goad 1995, 
Ph.D. thesis); we elaborate in §3.1. For simplicity, we ignore 
the effects of cloud-cloud shadowing in assuming an effec- 
tive bowl surface, but do account for self-obscuration of the 
bowl by the outer rim (ie. lines of sight which pass through 
the obscuring dusty torus) which for this geometry occurs 
at inclinations i > 45 degrees for a > 2 (the bowl geometry 
power law index). We thus propose the bulk of the broad 
emission lines to form in gas clouds lying along an effec- 
tive surface which spans between the accretion disc at small 
scale height and the TOR spanning a range in scale heights 
(see also Czerny & Hryniewicz 2011). The dimly illuminated 
gas clouds beyond the broad emission line region bowl sur- 
face are then likely dusty. The velocity field is described 
by equation 4 (discussed further in §3.3), and importantly 
we include the effects of Transverse Doppler Shift (hereafter 
TDS), Gravitational Redshift (hereafter GR), and turbu- 
lence on the emergent line profile (see A3. 21 and ^3.3)) . 

Initially, we parameterise the radial surface line emis- 
sivity distribution F(r) as a simple power- law in radius 
(F(r) oc r 7 ), with power- law index 7 = — 1 which is a fair 
approximation to the expected radial emissivity distribu- 
tion derived from photoionisation calculations for several of 
the commonly observed UV and optical emission- lines (e.g., 
figure [9]). We further assume that the line-emitting gas re- 
sponds linearly to variations in the incident ionising photon 
flux $h (i.e., locally, the marginal response of the line to 
continuum variations remains constant with time) . The spa- 
tially dependent line responsivity 77 (r), is here defined as in 
Korista & Goad (2004), i.e., 



AlogF(r)^ 
A log $ H 



(3) 



where F(r)u ne is the radial surface line emissivity, and <J>h 
is the incident hydrogen ionising continuum flux. For small 
amplitude continuum variations, this definition of responsiv- 
ity converges to that given in Goad et al. 1993. Since we 
observe a spread in delays for lines of different species (see 
e.g., Clavel et al. 1991), in any given system, we focus here 
on bowl-shaped geometries parameterised by a bowl-shape 
power- law index a in the range 1 ^ a ^ 2.5, noting that in 
geometries with a < 1, gas at larger radii will not have a di- 
rect view of the central continuum source, while gas at small 
radii on the near side is obscured from sight even at rela- 
tively small viewing angles, a consequence of the steepening 
of the bowl sides as R decreases (increasing dH/dR x ). 

In figures 3-6, we show the spatial distribution of 
BLR gas (as viewed from the azimuthal direction, R 2 ), 
the emissivity- weighted 2-d and 1-d response functions and 
steady-state line-profiles for bowl-shape power-law index a 
in the range 1.0 ^ a ^ 2.5 (with the bowl geometry normal- 
isation constant f5 chosen to match the requirement that the 
maximum time-delay at the outer radius when viewed face- 
on is 50 days) , and the fiducial geometry emission line profile 
as a function of the line of sight observer inclination, i =5, 
10, 30, and 45 degrees. At large line-of-sight inclinations, 
the surface of the bowl becomes increasingly self-obscured 
since such sight lines must first pass through the surround- 
ing dusty torus (the inclination at which self-obscuration 
first occurs decreases with increasing a). 

As expected the form of the 2-d response functions are 
broadly similar to those obtained for geometrically thin discs 
(the superposition of elliptical response functions) with a 
few notable differences. Firstly, as a increases, the fraction 
of surface area that lies at larger radii increases, so that the 
BLR response becomes increasingly weighted toward larger 
time-delays and thus lower line of sight velocities (e.g., fig- 
ure 4 left-right). Furthermore, since rings at larger radii, are 
elevated with respect to those closer in, there is a smaller 
offset in time-delay between the centre of each ellipse than 
would be expected for a geometrically thin disc. This eleva- 
tion of material out of the plane of the disc at larger radii, 
results in bowl-shaped geometries displaying enhanced re- 
sponse at small time-delays when compared to thin disc ge- 
ometries with similar Ri n , Rout- We note that for fixed R ou t, 
the surface area of the bowl decreases with increasing scale 
height if, and increases with increasing a. 

The increased weighting toward larger radii with in- 
creasing a is most evident in the 1-d response functions. 
For a = 1, the 1-d response function is essentially flat with 
increasing time-delay, falling after a time-delay R ou t(l — 
cos(0 — z)), corresponding to the light- crossing time of that 
part of the BLR which lies closest to the line of sight. At 
large a, the increased weighting of the emissivity to larger 
radii (an area effect) causes the 1-d response function to 



4 Here we chose n = 1.0 noting that for n(r,t) = constant, the 
size of n effects only the amplitude of the line response. In reality, 
77(7*, t) = constant is a poor approximation for most lines (see e.g., 
figure [9] of this work, and Goad et al. 1993, and Korista &; Goad 
2000, 2004). 
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Figure 3. The spatial distribution of BLR clouds as viewed from 
the polar axis (R z direction). Individual panels show the expected 
cloud distribution for bowl-shape power-law index a in the range 
1.0 ^ a ^ 2.5, and inclination i=5, 10,30 and 45 degrees, a = 1 
represents a special case wherein individual clouds populate the 
surface of a cone. For our model, no shadowing is assumed to 
take place. The grey-scale represents the observed intensity (here 
normalised to the peak), for our adopted power-law emissivity 
distribution F(r) oc r 7 , with power-law index 7 = — 1. 



peak at R ou t(l — cos(0 — moving towards smaller delays 
as i increases. In the lower right panel of figures 4 and 5 
the occultation of gas close to the line of sight results in the 
appearance of a secondary peak in the response function at 
small time-delays. 

In all cases, the emission-line profiles are double-peaked, 
and broadly similar to those found for thin disc geometries. 
The locations of the peaks are as for thin disc geometries 
determined by the velocity at the outer radius (as given by 
equation 2), and the inclination i (ie. zbvsini). Obscuration 
by the surrounding dusty torus at large line of sight incli- 
nations can be identified in the bottom right hand panels of 
figure 6, by the far larger difference in height between the 
line peaks and line centre and the U-shaped appearance of 
the line core, resulting from obscuration of low- velocity gas 
at large BLR radii and lying close to the line of sight. 

While the parameter space to be explored is clearly 
large, figures 4-6 indicate that the emissivity-weighted 2- 
d and 1-d response functions and steady-state emission-line 
profiles are broadly similar for fixed observer inclination. 
Thus for the remainder of this work, we explore the obser- 
vational characteristics of a single reference model, which 
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Figure 4. The corresponding 2-d response functions for figure [3] 



we refer to as our fiducial BLR geometry, parameterised by 
bowl-shape power-law index a = 2 (corresponding to a nor- 
malisation constant /3 — 1/150, for R and H measured in 
It- days). 



3.1 The role of emission line anisotropy and bowl 
scale height 

For isotropically emitting gas orbiting in a geometrically 
thin disc, the centroid of the 2-d response function and 
the emission-line profile 'shape' as quantified by the ratio 
fwhm/ai in the absence of transverse Doppler shift (TDS) 
and gravitational redshift (GR) are independent of inclina- 
tion. Furthermore, for our adopted line radiation pattern 
for individual clouds (which roughly speaking approximates 
a radiation pattern similar to the phases of the moon, e.g. 
O'Brien et al. 1994), increased anisotropy increases the cen- 
troid of the 2-d response function as i increases, but the 
profile shape remains unaltered (reduced emission on the 
side nearest to the observer due to increased anisotropy in 
the line is replaced by enhanced emission on the far-side at 
the same velocity - front back symmetry). 

By contrast for bowl-shaped geometries, the centroid of 
the 1-d response function increases with increasing observer 
line of sight inclination even for fully isotropically-emitting 
clouds (since rings on the bowl are pivoted about an axis 
running through the centre of the ring and the centre of 
the bowl). In addition, increased emission line anisotropy 
(see also Ferland et al. 1992) reduces the centroid of the 1- 
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d response function at small inclinations (by reducing the 
contribution of gas at high elevations relative to that at 
low elevations), while increasing the centroid at larger in- 
clinations due to the reduced contribution of gas closer to 
the observers line of sight. However, as for a thin disc ge- 
ometry, front-back symmetry of our adopted line radiation 
pattern, yields line profile shapes which are independent of 
inclination in the absence of GR, TDS effects. Notably, for 
bowl-shaped BLR geometries, the centroid of the response 
function is a stronger function of inclination when compared 
with that for a standard thin disc geometry. 

In practice we are mainly interested in geometries for 
which the line of sight inclination i allows an un-obscured 
view of the bowl surface (0 ^ i < <fi), that is, over the edge 
of the TOR. The BLR will be obscured on the near-side 
(closest to the line of sight) at larger inclinations either by 
self-occultation by BLR gas lying at larger radial distances, 
or by the surrounding TOR unless the line of sight to the 
central continuum source covering fractions of these compo- 
nents are themselves low[f|. Indeed, low line of sight inclina- 
tions are suggested by a number of studies. X-ray studies 
of Seyfert 1 galaxies suggest BLR inclinations of i < 30 de- 
grees (Nandra et al. 1997, Nandra et al. 1999, Tanaka et al. 
1995). Disc model fits to the double-peaked emission- line 
objects favour inclination angles of between 18-36 degrees 
(Eracleous et al. 1994). Furthermore, if the measured open- 
ing angles of ionisation cones observed in type 2 AGN are 
representative of AGN in general, then type 1 objects must 
be viewed at angles of i < 35-60 degrees (Wilson and Tsve- 
tanov 1994). 




4208642086420 
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Figure 5. The corresponding Id response functions for figure [3] 



3.2 Transverse Doppler shift and Gravitational 
Redshift 

At the inner radius chosen for our fiducial BLR geometry 
(200 R g ), the effects of transverse Doppler shift (hereafter, 
TDS) and gravitational redshift (hereafter, GR) have a sig- 
nificant effect on the emergent line profile, at low line of 
sight inclinations for the majority of lines, and in particular 
for those lines which form preferentially at small BLR radii. 
Using the standard formulism for each of these effects, we 
compare in ngureOthe 2-d response functions and emission- 
line profiles for our fiducial BLR geometry for inclinations 
i=5 and i=10 degrees, with (right panels) and without (left 
panels) the combined effects of TDS and GR (the 1-d re- 
sponse function is not shown here as it remains unaltered 
by the inclusion of these effects). 

TDS and GR act to enhance the red wing response at 
small time delays (at the expense of the blue- wing response), 
thereby creating a strong red-ward asymmetry in the line 
profile. This effect is largest for small line of sight inclina- 
tions, since for our model, in the absence of turbulence, the 
line of sight velocity is zero at i = 0, if the effects of TDS 



5 The dusty TOR might well extend to substantially greater scale 
heights in this geometry than that demarcated by the location 
of the grain sublimation radius (e.g., see figure 1). In this case 
obscuration of the central regions can occur at smaller observer 
inclination angles. It is also possible that such sight lines are only 
partially obscured by dusty clouds, as in Nenkova et al. 2008. 



and GR are not included. GR effects will be most impor- 
tant for lines which form at small BLR radii (deep within 
the gravitational potential). Similarly for our chosen velocity 
field effects due to TDS will be stronger in lines which form 
at small BLR radii (e.g. compare figure [8] panels 5 and 6), 
where the velocities are more extreme. As shown in figure [3 
the strong red-blue asymmetry introduced into the line pro- 
file is significant, as is the enhanced red wing response and 
diminished blue-wing response in the 2-d response function 
at small time-delays. At an inclination of 10 degrees, the 
importance of these effects diminishes (figure [3- right pan- 
els, cf. dashed line), and at 30 degrees are barely detectable 
(fwhm/ai ~ constant for i > 20 degrees). A key result of 
this work, is that the effects of GR and TDS acting alone 
can provide substantial width to the emission-line (~ sev- 
eral hundred km s _1 ), even for flattened, thin- or thick disc 
geometries, with pure planar Keplerian motion, viewed at 
low line-of-sight inclination. 

6 We do not consider here the effects of internal line broadening 
by electron scattering within the BLR clouds. For the physical 
conditions extant within the BLR gas, electron scattering optical 
depths of up to r e ^0.1 are expected (for cloud ionisation param- 
eters U > 0.1, and column densities Nu > 10 23 cm -2 ; see Laor 
2006). Electron scattering optical depths of r e < 0.5 can readily 
account for broadening of up to a few hundred km s _1 in the 
scattered fraction of the line photons emerging from such a cloud 
(perhaps contributing to reducing the number of required clouds 
for smooth broad emission line profile wings; Arav et al. 1997, 
1998). However, unless the electron scattering optical depths are 
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Figure 6. The corresponding line profiles for figure [3] 

In the absence of these effects line profile shapes as 
quantified by fwhm/ai are independent of inclination in 
thin disc geometries. However, as we show in §6.1, in the 
context of flattened BLR geometries, when present GR and 
TDS introduce a strong inclination dependence to the line 
profile shape at low line of sight inclinations which matches 
both qualitatively and quantitatively the observed correla- 
tion between fwhm/ai and ai among AGN reported in the 
literature (e.g. figure 3 of Collin et al. 2006). 

While Mannucci, Salvati and Stanga (1992) considered 
relativistic effects when modelling the 1-d response and 
emission-line profiles of disc- and nest-shaped BLR geome- 
tries, the effects of TDS and GR on the form of the 2-d 
response function have largely been overlookecQ However, 
a prominent red-wing response is a key observational feature 
of the recently recovered 2-d response function for the opti- 
cal emission- lines in the Seyfert 1 galaxy Arp 151 (Bentz et 
al. 2010b, their figure 4), though the response on the short- 
large over most of the BLR, its effect on the emission line profiles, 
even at small inclinations, should be minor. 

7 Eracleous &; Halpern (1994, 2003), considered relativistic ef- 
fects when modelling the double-peaked broad optical emission- 
line profiles of a sample of radio-loud AGN. Corbin (1997) con- 
sidered these effects on the shapes of the line profiles, but not 
on the reverberation line response functions. Kollatschny (2003b) 
argued that the redshifted component of the rms profile in the 
Seyfert 1 galaxy Mkn 110 was a result of gravitational redshift. 
From this he deduced a black hole mass of 1.4 x 10 8 Mq and an 
inclination of i=19 degrees for this object. 




-2000 2000-2000 2000 -2500 2500 -2500 2500 
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Figure 7. Upper panels - 2-d response functions for our fidu- 
cial BLR geometry, with (right-hand panels) and without (left- 
hand panels) the combined effects of Transverse Doppler Shift 
and Gravitational Redshift. Lower panels - the corresponding 1- 
d line profiles. The dashed line in the right-hand side of panel 
4 shows the same model without these effects (appearing in the 
left-hand side of panel 4). 

est timescales remain temporally unresolved. While these 
authors used the recovered response function to rule out a 
number of simple models for the BLR in this object, the 
effects of TDS and GR were not considered. If instead the 
enhanced red-wing response is a direct consequence of the 
combined effects of TDS and GR, then potentially Arp 151 
may be a system with a low line-of sight inclination, or a sys- 
tem with significant emission at small BLR radii, for which 
the effects of TDS and GR are both larger. 



3.3 Macro-turbulent cloud motion 

For material to accrete, significant angular momentum must 
be removed via viscous dissipation outward through the disc. 
However, this process does not require the disc to be geo- 
metrically thin. Indeed, the relative ionising continuum and 
broad emission-lines fluxes require covering fractions for the 
BLR gas of at least 10%, and often substantially larger. 
This means that even for thin disc BLR geometries, the 
disc must have significant scale height in order to intercept 
sufficient ionising continuum necessary to explain the ob- 
served line strengths. Furthermore, moderately flared discs 
are predicted by standard accretion disc models (Shakura & 
Sunyaev 1973), and turbulence remains the simplest mecha- 
nism for supplying the internal pressure necessary to support 
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the discs' vertical extent. While the predicted scale heights 
for thin-disc geometries are rather modest (H/R ~ 0.1), in 
our model the BLR is not physically attached to the disc, 
though indeed it may be the originator of the BLR clouds. 
Instead we model the BLR as an ensemble of line-emitting 
clouds, and invoke macroscopic turbulent cloud motion as 
the mechanism which provides the BLR geometry with the 
necessary scale-height to intercept a substantial fraction of 
the ionising continuum. Recent dynamical models of the con- 
tinuum and broad emission-line variations in the Seyfert 1 
galaxy Mrk 50 (Pancoast et al. 2011), suggest flared discs 
with significant scale height (~ 0.5 at the outer radius, for 
their assumed opening angle of 25±10 degrees) as a plausible 
geometry for this system. We also note that the surround- 
ing TOR, which is the likely reservoir of material feeding 
the black-hole, itself must have substantial scale-height in 
order to obscure the BLR in type 2 objects (in the general 
orientation-dependent unified picture for AGN). 

While GR and TDS provide substantial line width in 
thin disc geometries even for face-on objects, the line width 
is still largely determined by the black hole mass and source 
inclination. Thus for NLSyls, which are thought to be low 
mass high accretion rate systems observed at low inclina- 
tions, a significant turbulent component may still be re- 
quired to ensure sufficient width in the line0 

In order to implement isotropic turbulence within the 
context of our model we add in quadrature a randomised 
velocity component Vt ur b which increases the local Keplerian 
velocity v^ ep (equation 2) according to 



v 



2 

: ^kep 



i+[^y 



(4) 



and such that |iWb| increases linearly with scale height H, 
i.e., | vturb | = frturb Vkep, where b tU rb is a constant of 

order unity. This is similar to the implementation of turbu- 
lence with BLR scale height suggested in Collin et al. (2006). 
In practice, we draw at random a randomised (in direction) 
velocity vector from a Gaussian distribution of the appro- 
priate width. Using this formulism, the turbulence is zero 
both when frturb = (turbulence switched off) and when 
H — (zero scale height). Thus in the context of bowl- 
shaped BLR geometries, the contribution that turbulence 
makes to the 2-d response function and line profile depends 
upon the bowl-shape, and on the radial surface line emissiv- 
ity distribution which together determine the scale height 
at which a given line forms. Thus, in our model, the effects 
of turbulence are minimised for lines formed near the base 
of the bowl (small scale height). For non-zero H and suffi- 
ciently large 6 tU rb, the randomised turbulent velocity com- 
ponent dominates over the planar Keplerian motion. 

In figure [8] we illustrate the effects of turbulence on 
the 2-d response function and emission-line profile for our 
fiducial BLR geometry assuming a power-law emissivity 
distribution with power- law index 7 = —1, and including 
the effects of TDS and GR, for an inclination i=5 degrees 
and turbulence parameter 6 t urb = 1 (left-hand panels), and 



8 We speculate that high Eddington rate sources may have TORs 
with larger than typical covering fractions, necessitating smaller 
observed viewing angles i for un-obscured lines of sight. 



frturb = 2 (middle 2 panels). Turbulence has a number of im- 
portant attributes. First, the lines are broader (as expected), 
and at low inclinations the prominent shoulders/horns in- 
dicative of flattened BLR geometries have disappeared even 
for moderate turbulence. Secondly, gas at large radii now 
makes a significant contribution to the line wings, enhanc- 
ing the line wings relative to the line core. The effect of 
turbulence on line shape depends upon the size of b tU rb and 
on the line formation radius. For example, as expressed in 
terms of the typically measured parameters used to quan- 
tify line profile shapes, the ratio of the line full- width at half 
maximum (hereafter fwhm) to the line dispersion (hereafter 
ai, see e.g., Collin et al. 2006), then for b tU rb = 0, 7 = — 1 
(dotted red line, lower left panel), fwhm/cri = lA, increas- 
ing to 1.5 for frturb = 1 (solid black line). For larger 6turb, 
fwhm/ cri = l.l (lower middle panel, solid black line). While 
in panel 3, where we have now increased the radial power- law 
emissivity to 7 = — 2, fwhm/ai increases from 1.0 without 
turbulence (dotted red line) to 1.4 with 6 tU rb = 2 (solid black 
line) . 

We note that the observed lack of response of the broad 
emission-line wings on short time-scales in some objects has 
traditionally been explained by the presence of a low respon- 
sivity component in the inner BLR (see e.g., Goad et al. 
1993; Sparke 1993; Korista & Goad 2004). Here, a similar 
effect arises for a different reason, a significant fraction of 
the broad emission-line wings, as observed at low inclina- 
tions, are now due to gas lying at large BLR radii which 
responds on longer time scales. If turbulence dominates the 
contribution to the emission-line response, then the core and 
wings of the line will respond on similar timescales. 

One potentially significant result of this work is that at 
low inclinations, and in the presence of turbulence our model 
line profile displays a striking resemblance to a Lorentzian 
function. By way of illustration, we model the line profile 
with a 3-parameter Lorentzian of the form 



(a 



■ xo) 



+ r 2 



(5) 



where Io is the peak intensity, xo is the median velocity, 
and r represents the half-width at half maximum of the 
line. Our best fit model shown in the lower middle panel 
of figure [8] (dashed blue line) with fit parameters, I x =o = 
1.025, T = 623.5 km s~\ and x = +82.2 km s" 1 provides 
a remarkably good fit to the line profile, failing only in the 
extreme line wings which would in any case be difficult to 
measure in emission- line spectra due to their low contrast. 
The appearance of narrow cores and extended broad wings 
as observed in some type 1 Seyferts (in particular NLSyls 
or "Pop A" objects; see Collin et al. 2006) may indicate the 
presence of significant turbulence and low inclination (< 20 
degrees) in these objects. While many studies have modelled 
the Balmer and Mgn emission line profiles with Lorentzians 
(e.g, Zamfir et al. 2010, their "Pop A" type 1 AGN), to 
our knowledge no explanation of why such profiles should 
manifest themselves has ever been offered, heretofore. 

For our model, turbulence is more significant for lines 
formed at larger BLR radii (e.g., H/3, Mg 11), and hence 
larger scale heights. This is illustrated in figure [8] - right- 
hand panels, where we show the effect of steepening the 
power- law emissivity distribution to a slope of —2. The bulk 
of the line is now formed at small BLR radii, and thus small 
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Figure 8. Upper panels - 2-d response functions for our fidu- 
cial BLR geometry with a line-of-sight inclination 5 degrees and 
turbulence parameter 6 tur b = 1, (left panel), 6 tur b = 2 (mid- 
dle panel). In panel 3 we show 6turb — 2 for a power-law emis- 
sivity distribution with 7 = —2. Lines formed at smaller scale 
height (with steeper emissivity distributions) are less affected by 
turbulence. Also shown (lower panels) are the steady-state line 
profiles (solid black lines), together with the expected profile for 
no turbulence, 6 tur b = (red dotted line). In the middle lower 
panel we also show our best fit Lorentzian profile (blue dashed 
line) with fit parameters, I = 1.025 , T = 623.5 km s -1 , and 
x = +82.2 km s _1 . 



domising the velocity field) . If turbulence is significant these 
lines' emission-line profiles will show enhanced wings rela- 
tive to the line core. For low inclination systems, turbulence 
produces Lorentzian profiles, especially in lines forming at 
large radii (e.g., Balmer and Mg 11) and so substantial scale 
height, similar to those seen in many NLSyl (or "Pop A") 
spectra. At higher inclinations, the Keplerian velocity domi- 
nates over the turbulent component, although the latter still 
affects especially the cores of the emission-line profiles (see 
Figures Al, A2). In light of these findings, we suggest that 
the previous association of an enhanced red-wing response 
with in-falling gas or a Keplerian disc + hotspot (for exam- 
ple, as suggested for the velocity field in the nearby NLSyl 
Arp 151 (Bentz et al. 2010) and indeed spherical or disc- 
wind models of the BLR (Konigl and Kartje 1994; Chiang 
and Murray 1996; Murray and Chiang 1997), may also in 
part be attributed to the presence of mildly relativistic gas 
in the inner BLR, as we have demonstrated here, as was 
previously suggested to account for the enhanced red-wing 
response observed in broad B./3 for the Syl galaxy Mrk 110 
(Kollatschny 2003b). 

We note here that the recovered 2-d response function 
and emission- line profile for R/3 in the NLSyls NGC 4051 
(Denney et al. 2011, their figure 4), Arp 151 (Bentz et al. 
2010, their figures 2,3 and 4), and in the 2-d CCF of Mkn 110 
(Kollatschny 2003, figure 7), bares an intriguing resemblance 
to that predicted for a bowl shaped geometry viewed close to 
face-on with a significant scale height dependent turbulent 
component (e.g., compare their figures with the left hand 
panel of Figures A1-A3 in Appendix A, for i = 10 degrees). 
In Arp 151, an enhanced red-wing response is evident in all 
of the Balmer lines. 

Significantly, GR, TDS and turbulence can on their own 
provide substantial line-width even for purely transverse mo- 
tion and may explain the observed line widths in systems 
generally thought to be viewed close to face-on as well as the 
break in the fwhm/ai relation at small o~i observed among 
the AGN population. 



scale heights where the effect of turbulence is low. Conse- 
quently, though the turbulence parameter is large, a strong 
red-blue asymmetry, a result of TDS and GR whose effects 
are larger at small BLR radii, can still be discerned in the 
2-d response function and emission-line profile. Thus at low 
line of sight inclinations, we expect to observe strong differ- 
ences between the broad emission-line profiles of low- and 
high-ionisation lines. 

In summary, TDS, GR and turbulence play a signifi- 
cant role in determining the shape of the 2-d response func- 
tion and emission-line profile. Lines formed preferentially at 
small BLR radii will show an enhanced red-blue asymme- 
try in their 2-d response function and line profile due to the 
combined effects of TDS and GR (these effects will remain 
apparent provided the emission from the turbulent compo- 
nent does not dominate the contribution to the line profile) . 
By contrast lines formed at large BLR radii will display more 
symmetric 2-d response functions and line profiles (because 
the effects of GR and TDS are diminished at larger BLR 
radii, and turbulence at larger scale height acts to reduce 
the dependence of line profile shape on inclination by ran- 



4 SIMULATIONS 

4.1 Photoionisation calculations 

As in previous work, we adopt the most general model for the 
BLR gas, the locally optimally emitting cloud model of Bald- 
win et al. 1995 (see also Korista & Goad 2000, 2001, 2004). 
In brief, this model assumes that there exists a large popu- 
lation of BLR clouds with a broad range of gas density, gas 
column density, and ionising source distance. The observed 
broad emission-line spectrum arises by a process of natural 
selection, since for a given line only clouds whose properties 
span a relatively narrow range of physical conditions are ef- 
ficient in re-processing the incident ionising continuum into 
particular line emission (hence "locally optimally emitting" , 
or LOC). This model not only explains why the observed 
spectra of AGN are generally similar, but significantly, re- 
quires minimal fine-tuning. 

Using the photoionisation modelling code CLOUDY 
version c08.00 (Ferland et al. 1997, Ferland et al. 1998), 
we generated a grid of photoionisation models for simple 
slabs of gas (hereafter, clouds), each with a constant density 
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Figure 9. (a) Radial emissivity distributions for selected low and 
high ionisation lines. The red line indicates a power-law emissivity 
distribution with power-law index 7 = — 1. (b) Radial surface re- 
sponsivity distributions for selected low and high ionisation lines. 
The dashed vertical lines indicate the BLR inner and outer radius 
for our fiducial BLR geometry. 



and an un-obscured view of the ionising continuum source. 
In this simple approach contributions from the reprocessed 
continuum (the diffuse component from the cloud) and the 
effects of cloud shadowing are ignored. Moreover, we do not 
consider the complicating effects of micro-turbulence, which 
would tend to broaden the intrinsic line profile far above 
the thermal width and reduce the central line optical depth, 
thereby increasing the line escape probability isotropically 
(in this case) and hence altering the emergent spectrum 
(see e.g., Bottorff et al. 2000). This would also act to re- 
lax the requirement of large cloud numbers to explain the 
relative smoothness of the observed broad emission-line pro- 
files (Arav et al. 1997; 1998). Of the four emission lines con- 
sidered here (see Figure [9]), micro-turbulence would affect 
mainly H/3 and Mg II, effectively moving their luminosity- 
weighted radii to smaller values. 

The full grid spans 7 decades in the gas hydrogen 
number density-hydrogen ionising photon flux plane, 7 < 
log 10 n H (cm~ 3 )< 14, and 17 < log 10 $ H (cm -2 s _1 )< 24, 
stepped in 0.25 decade intervals in each dimension (see 
also Korista et al. 1997). Since the emitted spectrum is 
not all that sensitive to the cloud column density over 
the range 22 < log 10 Ah (cm -2 )< 24, all clouds are as- 
sumed to have a constant total hydrogen column density 
log 10 A^H(cm _2 )= 23. The adopted elemental gas abun- 



dances are solar. The adopted spectral energy distribution of 
the incident continuum is a modified form of the Mathews & 
Ferland (1987) generic AGN continuum^ For the adopted 
mean ionising continuum source luminosity of NGC 5548, 
log 10 Li on — 44.14 erg s _1 , a hydrogen ionising photon flux 
log 10 $ = 20.0 photon/s, corresponds to a cloud-ionising 
source distance of R = 15 (Ho/70 km/s) light-days. 

Utilising the standard LOC gas density distribution 
(weighting) functional <7(tih) oc n^ 1 , spanning the range 
in gas density 8 < log 10 nu (cm _3 )< 12, as described in 
Korista & Goad (2000; see also Baldwin et al. 1995; Bot- 
torff et al. 2002), we sum the emission along the density 
axis for each radius, producing a radial surface line emissiv- 
ity function (erg cm -2 s _1 ) for each of the lines considered. 
The radial surface line emissivity distributions are computed 
beyond the point at which grains will form, i.e., beyond 
the point where the grain temperature at the illuminated 
face of the cloud falls below the grain sublimation temper- 
ature. The radius at which this occurs (~ 100 lt-days for 
our adopted continuum and the hardiest of grains) marks 
the inner boundary of the dusty TOR. For simplicity, we 
will assume the TOR to be opaque to UV/optical photons, 
and thus the central regions are obscured along lines of sight 
that pass through the TOR. The radial surface line emissiv- 
ity distributions for C iv A1549, Mg 11 A2798, H/3 A4861 and 
He 11 A4686 whose emission line profiles are commonly used 
in measuring black hole masses, are shown in figure EJa) 
(solid black lines). For our fiducial model we assume a con- 
stant inward fraction for each of the lines (approximately a 
luminosity weighted- average) of 70% for C IV and He II, and 
80% for Hf3 and Mg II, similar to the values determined for 
these lines over a broad range of physical conditions (e.g. 
O'Brien et al. 1994; Goad 1995, Ph.D. thesis; Korista et al. 
1997a). For individual BLR clouds, we adopt a line radi- 
ation pattern which approximates the phases of the moon 
(see e.g., O'Brien et al. 1994). 

In figure [9{b) we show the corresponding radial surface 
responsivities rj(r) (Eq. 3) for each line. These are derived 
from the radial surface line emissivities F(r), by measuring 
locally the slope of this distribution and dividing by —2. 
Since the radial surface line emissivities are a function of 
the ionising photon flux and our choice of weighting in the 
3>h, plane, so too are the radial surface line responsiv- 
ities. However, some general trends do exist. For example, 
the responsivities of H/3 and Mg 11 are anti-correlated with 
the incident ionising photon flux, being stronger at low inci- 
dent fluxes (in the outer BLR, where the optical depths are 
lower). This behaviour is also expected for He 11 within the 
inner BLR (see figure[9]). The C IV responsivity tends to align 
along lines of constant ionisation parameter (anti-correlating 
with U — $h/^hc), but also gently declines for larger inci- 
dent fluxes at constant U for gas densities nu ^ 10 10 cm -3 
(see also Korista & Goad 2004, their figure 4). 



9 The effects of a polar-angle dependent ionising continuum 
shape on the detailed emission line spectrum are explored in Ap- 
pendix [C] 

10 This standard LOC distribution (weighting) function in hydro- 
gen gas density is similar to that found by Krause et al. (2012), for 
BLR clouds confined by a magnetic field undergoing filamentary 
fragmentation. 
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4.2 The driving continuum light-curve 

A model that describes the BLR gas distribution, must be 
able to reproduce the measured broad emission-line intensi- 
ties and importantly their relative response to variations in 
the ionising continuum. While a single spectrum provides a 
"snapshot" of the physical conditions within the BLR, it is 
not a "true" representation of the current state of the line- 
emitting gas. This arises because the BLR is physically ex- 
tended and therefore a single spectrum comprises contribu- 
tions from gas in response to prior continuum states (due to 
the finite light-travel time across the BLR) up to and includ- 
ing the continuum state at the current epoch. Consequently, 
in order to compare the general characteristics of our fiducial 
bowl-shaped BLR model (see Table 1) with observations, we 
must first drive it with a model continuum light-curve whose 
variability power is representative of that determined for the 
source in question, taking into account both sampling effects 
and the finite duration of the light-curve (data sampling can 
be so sparse that the highest frequencies present in the un- 
derlying PSD are only poorly sampled by the observations). 
Our aim here is not to match the detailed response of indi- 
vidual lines, nor to fit the emission-line spectrum at a single 
epoch (for the reasons outlined above), rather here we aim 
to explore the gross properties of bowl-shaped geometries 
by determining the time-averaged spectrum and the aver- 
age emission-line response timescale for the strongest UV 
and optical recombination lines for comparison with obser- 
vations. 

Here we adopt the approach of Kelly et al. (2009), 
and model the driving continuum light-curve as a damped 
random walk (see Kelly et al. 2009, and Macleod et al. 
2010, for details), which has been shown to be an ex- 
tremely good match to the variability of quasars (Kozlowski 
et al. 2010). We set the characteristic continuum variabil- 
ity timescale r c har to that measured by Collier and Peterson 
2001, Tchar = 40 days for NGC 5548, with a total duration 
of 450 days and daily sampling. The duration of our light 
curve is chosen to approximately match that of a typical 
monitoring season noting that the first r m ax points must 
be discarded since at prior times, the outer BLR has yet to 
respond to the continuum variations. 

For each input continuum light-curve we compute the 
corresponding velocity resolved emission-line light-curve 
L(v,t), assuming that locally, the emission-line gas responds 
linearly to continuum variations, i.e. a locally-linear re- 
sponse approximation (see Goad et al. 1993), so that while 
77 (r) remains time-independent, we now take into account 
the radial variation in line responsivity r/(r). As shown in 
Figure [9jb), and by Goad et al. 1993, the emission- line 
line responsivity shows a strong radial dependence. How- 
ever, a locally linear response has been shown to be a rea- 
sonable approximation provided that the amplitude of the 
continuum variations are small (Goad, 1995, Ph.D. thesis, 
O'Brien et al. 1995). For larger continuum variations the lo- 
cally linear response approximation is no longer valid and 
a full treatment of the time-dependent non-linear effects is 
required 

11 We do not consider changes in the emission-line response due 
to variations in ionising continuum shape nor variations in the 
inner and outer BLR boundaries. In any case, it is likely that local 



For each combination of continuum-emission-line light- 
curve we compute the peak and the centroid of their cross- 
correlation function (CCF, see e.g., Gaskell & Peterson 
1987), for the latter adopting a threshold of 0.6 for the 
centroid calculation, as well as four different measures of 
the velocity field, the fwhm and dispersion of the mean 
and root-mean square (hereafter rms) profiles, quantities 
that are commonly used by observers in the estimation of 
black hole masses using the reverberation mapping tech- 
nique and/or scaling relations derived therein. To ensure 
that the full range in light-curve behaviour is sampled we 
repeat this process 1000 times and thereafter compute prob- 
ability distribution functions for each of the measured quan- 
tities. In the following section we present the results of our 
simulations and place them within the context of results re- 
ported in the literature. 



5 RESULTS 

Before exploring the effects of reverberation within the phys- 
ically extended BLR on measurements of the continuum- 
emission-line delay and emission-line velocity dispersion, pa- 
rameters typically used in virial mass estimates of the cen- 
tral black hole, we start by investigating an idealised case in 
which the predicted delay and velocity dispersion have been 
measured from the instantaneous 1-d responsivity- weighted 
response function and variable line profile. Such a situation 
would arise, for example, if the BLR were illuminated by 
a delta-function pulse rise in the ionising continuum flux. 
In so doing, we aim to reveal the relationship between the 
physical properties of our model (ie. our assumed geometry, 
emissivity, responsivity, velocity field, turbulence, and incli- 
nation) and the measured quantities used in black hole mass 
determinations. 

The black hole mass estimates are here based on a virial 
product formed from the mean square dispersion, cr 2 , of 
the variable line profile and the centroid, r ce nt, of the 1-d 
responsivity- weighted response function, such that. 



where Rbw — cr ce nt, and 5v 2 = af. Since the black hole mass 
is an input to our model, it is straightforward to determine 
the virial scale factor / required to force agreement between 
the measured black hole mass and the input valuJ* 2 !. 

We have calculated virial masses, virial scale factors /, 
and line shapes for two pairs of emission-lines, H/3-Mg II, 
and C iv-He 11, chosen because they are representative of 
lines typically measured in ground- and space-based AGN 
monitoring campaigns, and importantly because they span a 

variations in n(r) due to changes in ^(H) will be the dominant 
non-linear effect for the majority of lines. A thorough treatment 
of these non-linear effects in the context of bowl-shaped BLR 
geometries will be explored elsewhere. 

12 For our fiducial BLR geometry, a = 2, f3 = 1/150, the denom- 
inator in equation 2 is 53% larger ((4/3) 3 / 2 ) than for a standard 
Keplerian velocity field at the outer radius. Hence in the absence 
of turbulence, the velocities are lower than in the standard Ke- 
plerian model, and thus the masses derived using equation [6] are 
also lower, even before the effects of inclination are taken into 
account. 
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broad range in their characteristic line formation radii, and 
thus will probe the largest range in variability behaviour 
for our fiducial BLR geometry. We use the calculated radial 
surface line emissivity distributions of §4.1 (see figure 9), and 
adopt a locally-linear response approximation. We include 
the effects of anisotropic line emission, as described near 
the end of §4.1 

5.1 The importance of turbulence, line emissivity, 
and inclination on Mbh determinations 

To assess the impact of turbulence on the derived black hole 
masses and virial scale factor /, we have estimated Mbh for 
our fiducial BLR geometry for each of the four lines, for line- 
of-sight inclinations spanning the range 2-50 degrees and 
using 3 values for the turbulence parameter, b tUT h = (no 
turbulence), frturb = 1 (moderate turbulence) and 6turb = 2. 
In figure [10] we show the virial mass estimates using equa- 
tion [6] for each of the four emission lines (upper 4 panels), 
as well as the virial scale factor / (middle 4 panels), re- 
quired to reproduce the input black hole mass, as a func- 
tion of line-of-sight inclination i. Individual colours indicate 
the degree of turbulence, frturb = (black line), frturb = 1 
(red line) and frturb = 2 (green line). We find that in all 
cases, the estimated black hole mass is lower than the input 
mass, and is a strong function of inclination, being larger 
at larger line-of-sight inclinations. The middle 4 panels of 
figure [TOl indicate that for a given inclination, and in the ab- 
sence of reverberation effects, each line requires a different 
virial scale factor to recover the input mass. The discrepancy 
between the emission-lines is largest for zero turbulence and 
low line-of-sight inclinations. Additionally, black hole mass 
estimates are systematically lower and thus /-factors sys- 
tematically higher with turbulence switched off (solid black 
line) . The discrepancy between the virial mass estimate and 
the input mass can be larger than 2 dex for near face-on 
systems with turbulence switched off. For lines which form 
at large BLR radii, corresponding to a larger scale height in 
our model, e.g. H/3, Mg II, turbulence significantly enhances 
the black hole mass estimates at low line of sight inclina- 
tions (solid green line), reducing the discrepancy between 
the virial mass estimate and the input mass by more than 1 
dex. This arises because for our chosen velocity field, in the 
absence of turbulence, the line width at small line-of-sight 
inclinations is a result of the combined effects of TDS and 
GR only. Turbulence in our model acts to increase the lo- 
cal velocity field at large scale heights thus broadening the 
emission-line and thereby increasing the virial mass estimate 
and consequently lowering the virial scale factor. Since in our 
model turbulence increases with increasing scale height this 
effect is more pronounced in lines which form at large BLR 
radii (e.g. H/3, Mg n). Furthermore, turbulence randomises 
the direction of the velocity field and thereby acts to reduce 
the otherwise strong dependence of the emission-line width 
on inclination, normally found for planar Keplerian motion. 

In the lower panel of figure [10] we show the emission- 
line shape (fwhm/ai) as a function of line fwhm (or equiv- 
alently line-of-sight inclination i). This reveals a strong de- 
pendence of emission-line shape on line fwhm, such that 
broader lines have more Gaussian profiles, while narrower 
lines are more Lorentzian in shape. In general turbulence 
acts to soften the otherwise strong dependence of line shape 
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Figure 10. Panels 1-4 - Mbh as a function of line of sight incli- 
nation i for each line for our fiducial BLR geometry and 3 values 
of the turbulence parameter, 6 tur b = (black line), 6 tur b = 1 (red 
line) and 6 tur b = 2 (green line). Mass estimates are here based on 
measurements of the mean-square dispersion of the steady-state 
(mean) line profile and centroid of the 1-d responsivity-weighted 
response function. Panels 5-8 - the corresponding virial scale fac- 
tors, / also as a function of inclination i. Panels 8-12 - depen- 
dence of line shape fwhm/ai as a function of line fwhm (for the 
same 3 values of &turt>)- 



on inclination. However at low line-of-sight inclinations, tur- 
bulence dominates over the planar Keplerian motion moving 
low velocity gas from the line core to the line wings, result- 
ing in profiles with narrower cores and extended line wings 
(fwhm/ai < 1). Once again, this effect is most pronounced 
in lines which form at large BLR radii (e.g. H/3, Mg II). Fig- 
ure [10] indicates that for lines formed at small scale heights 
(C IV and He il), the effect of turbulence on the derived Mbh 
and /-factors is significantly smaller. 

Onken et al. (2004) showed that nearby AGN for which 
both reverberation based and velocity dispersion based mass 
estimates are available, follow the same Mbh — cr relation as 
quiescent galaxies. From this they derived a statistical esti- 
mate of the virial scale factor for all AGN with reverberation 
mapping data, finding an average /-factor of / = 5.49 ± 1.65 
for H/3, see also Table 2 of Collin et al. (2006), here indicated 
by the horizontal solid and dashed lines (middle panel of fig- 
ure [TO}. Woo et al. (2010) provide an updated value for the 
virial scale factor /. Their estimate log/ = 0.72 zb 0.1 is 
based on matching reverberation masses of 24 AGN with 
those obtained from recent stellar velocity dispersion esti- 
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Figure 11. Panel 1 (upper left) - Mbh as function of inclination 
for a geometrically thin disc and power-law emissivity distribution 
with power-law index 7 = (blue lines/symbols), —1 (green) and 
—2 (red). Also shown (green dashed lines) is the effect of changing 
the line radiation pattern (for 7 = — 1) from 50% (isotropic)- 
100% (anisotropic). Panel 2 (lower left) - corresponding virial 
scale factors. Panel 3 (upper right) - corresponding line profile 
shapes as a function of fwhm. Panel 4 (lower right) fwhm/ai as 
a function of power-law emissivity index 7. Differences in the ratio 
fwhm/ a 1 between lines may result from differences in their radial 
surface emissivity distributions. N.B. a geometrically thin disc has 
a very small scale height by construction, and thus the turbulent 
velocity is effectively zero even though we chose 6turb = 2. 



mates of the host galaxy, adopting the Mbh — cr relation of 
quiescent galaxies from Giiltelkin et al. (2009), and is con- 
sistent within the errors with that found by Onken et al. 
(2004). We note that for our fiducial BLR geometry, the 
inclusion of a scale height dependent turbulent component, 
increases the range in inclination returning viral scale fac- 
tors that are consistent within the errors with the average 
/-factor for H/3 reported by Collin et al (2006) and Woo et 
al. (2009). 

For completeness we show in the left hand panels of fig- 
ure [TT] the derived black hole mass and virial scale factor, 
/, as a function of inclination for a geometrically thin disc 
with the same inner and outer radius as our fiducial BLR 
geometry (1.14 and 100 light-days respectively), turbulence 
parameter b tU rh = 2, and using simplified power- law line 
emissivity distributions with power-law index 7 = (solid 
blue line), —1 (solid green line), and —2 (solid red line). We 
also show (dashed green lines) the effect of changing the line 
anisotropy from 50% (ie. isotropic, Fi nw d = \F to ti) to 100% 
(ie. anisotropic, Fi nw d = F to ti)- Line anisotropy increases 
the emissivity- weighted radius by a factor 1 + \ sin 2 i, (Goad 



1995, Ph.D. thesis, O'Brien et al. 1994), and thus for the lim- 
ited range of inclinations studied here results in a marginal 
increase in the virial scale factor for a geometrically thin 
disc. As for our bowl-shaped geometry the masses are sys- 
tematically underestimated (except at large inclination and 
for steep emissivity distributions - red line, upper left-panel) 
and the discrepancy between the calculated and input black 
hole mass decreases with increasing inclination. The thin 
disc model also shows a significant change in profile shape 
at small inclinations (figure ITU upper right panel). Since 
a disc by construction has little scale height, the contribu- 
tion of the turbulent component is minimal (even though 
frturb = 2). Thus for a given radial surface line emissiv- 
ity distribution, only GR and TDS can play a significant 
role in modifying the line shape for thin-disc geometries at 
small inclinations. We note that in the absence of these ef- 
fects, the emission-line shape is independent of inclination 
for our adopted emission-line radiation pattern. Of more sig- 
nificance is the strong dependence of the mass estimates on 
the radial surface line emissivity distribution, which indi- 
cates that in the absence of turbulence, lines which are pref- 
erentially formed at small BLR radii yield larger black hole 
masses, and thus smaller virial scale factors. A similar re- 
sult was found for our fiducial model (figure fTOl— solid black 
lines, cf. H/3-Mg 11 and C iv-He 11). We discuss this further 
in the following section. 

Finally, we note that unless turbulence dominates the 
velocity field, then for a bowl-shaped BLR geometry, we ex- 
pect to observe a strong f-fwhm dependence in all lines re- 
gardless of where they form (due to the strong dependence of 
fwhm on inclination), and broadly similar to that found for 
disc- like configurations (see Decarli et al. 2008, their figure 
6). Scatter in the relation between different lines may point 
to differences in their radial surface line emissivity distribu- 
tions and/or the presence of a significant turbulent velocity 
component. 



5.2 fwhm/ai - the role of emissivity 

As noted by Collin et al. (2006) the ratio of the zeroth and 
2nd order moments (mean and root mean square respec- 
tively) of the broad emission-line profile can be used to place 
constraints upon the BLR geometry and kinematics. Decarli 
et al. (2008) found that in a sample of 36 AGN with roughly 
equal numbers of radio-loud and radio-quiet objects, that 
the ratio fwhm/ai for H/3 is closer to that derived for an 
isotropic BLR geometry (wherein fwhm/ai = 2.35), while 
the smaller values (~ 1) measured for C IV are more sugges- 
tive of a flattened BLR geometry. Combined with the strong 
correlation between their virial scale factor fj* 3 ! and fwhm 
found for both lines (such a correlation is expected if line 
of sight velocity is a strong function of inclination) , and the 
absence of a correlation between the fwhm of H/3 and C IV, 
Decarli et al. (2008) argue that their results are consistent 
with C IV originating in a flattened BLR geometry, with H/3 
originating in a geometrically thick disc with a significant 



13 Decarli et al. (2008) use a different definition of the virial scale 
factor fd, ^blr — fd x fwhm, which using our definition implies 
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turbulent component. As in our model, the increased turbu- 
lence at large scale height reduces the dependence of the H/3 
line width on inclination and therefore could account for the 
reported absence of a correlation between the H/3 and C IV 
emission-line widths, the more isotropic appearance of the 
H/3 line profile, and if the turbulence is large enough, the 
larger width of H/3 relative to C I\0 

However, we caution here that the ratio fwhm/ai for 
fixed Rin , Rout has a strong dependence on the radial surface 
line emissivity distribution. For example, if we approximate 
the radial surface line emissivity distribution as a power-law 
in radius (F(r) oc r 7 ), then for a geometrically thin disc, the 
ratio fwhm/ai is a minimum for 7 = — 2, and increases for 
both smaller and larger 7 (see figure [11] lower right panel) . 
For 7 large and positive, the emission is weighted toward 
the outer edge only, yielding fwhm/ai « 2.8 as appropri- 
ate for a ring-like distribution. Similarly, for negative 7, the 
inner radius dominates the emission, while emission from 
rings at larger radii (and hence lower velocity), tend to fill 
in the dip at line core (between the horns) producing more 
rectangular- looking profiles (smaller fwhm/ai). In extreme 
cases (7 < —3), this can lead to fwhm/ai > 3. For our fidu- 
cial model, C IV has a steeper emissivity distribution than 
H/3 and thus a smaller fwhm/ai (cf. figure [TUl- lower panels 
1 and 3, and figure ITT1 - upper right panel), which may in 
part explain the differences found by Decarli et al.(2008). 
A similarly strong dependence of line shape on radial sur- 
face line emissivity distribution can be found for other BLR 
geometries (see also Robinson 1995a,b)E3 Thus while dif- 
ferences in the ratio of fwhm/ai between H/3-C IV may 
indicate differences in their scale height, we suggest that 
differences in their radial surface line emissivity distribu- 
tions likely also play a significant role. Evidence in support 
of this claim comes from the large range (factor of a few) in 
fwhm/ai displayed by the H/3 line in NGC 5548 during 13 
years of monitoring (Collin et al. 2006). Since in an individ- 
ual source neither the black hole mass nor its inclination can 
change appreciably on such a short timescale, this suggests 
that the observed variations in the Balmer line profile shape 
for NGC 5548, over the 13 years of observation, are due to 
gross changes in the radial surface line emissivity distribu- 
tion in response to large variations in the ionising photon 
flux within the physically extended BLR, though we cannot 
rule out dynamical changes on longer baselines (^several 
years for NGC 5548). 



5.3 The importance of Reverberation: Which 
measures of Rbw — Sv should be used? 

In order to assess the impact of continuum-emission-line re- 
verberation on the virial mass estimates and hence scale fac- 
tors /, we have driven our fiducial BLR geometry, with sim- 



14 We note that the number of objects for which simultaneous 
line width comparisons between H/3 and C iv have been made is 
small, and that claims of a non-correlation between their respec- 
tive line widths may be premature. 

15 In support of this claim, we show in Figure A2 the mean 
responsivity-weighted (black solid line) and emissivity-weighted 
(red dotted line) emission-line profiles, for both the low- and high- 
ionisation emission lines, for our fiducial model with turbulence 
parameter 6 tur b — 2, and a range of line of sight inclinations. 
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Figure 12. H/3 (hashed lines) and He II CCF centroid (left pan- 
els) and peak (right panels) distribution functions (resulting from 
1000 simulated continuum-emission-line light-curves) for a range 
of BLR inclinations. Also shown is the centroid of the 1-d re- 
sponse function for H/3 (vertical dashed line) and He II (vertical 
dotted line) as determined from the steady-state model. 



Table 1. Our fiducial BLR model. 
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ulated continuum light-curves whose variability characteris- 
tics have been designed to match those observed in the best 
studied AGN, the Seyfert 1 galaxy NGC 5548. Our fiducial 
BLR model (see Table 1 for details) utilises the computed ra- 
dial surface line emissivity distributions for each line (§4.1), 
our prescription for emission-line anisotropy and line radia- 
tion pattern (§4.1), and a turbulence parameter frturb = 2. 
All simulated emission-line light curves have been calculated 
assuming a locally-linear response approximation. For each 
of the four lines considered (H/3, Mg II, C IV and He 11), we 
compute their velocity resolved time-variable emission-line 
light-curve L(v,t), and their mean and root mean square 
(rms) emission-line profiles. 

Simulations were performed for a range of line of sight 
inclinations i=2, 5, 10,20,30,40, and 50 degrees, and repeated 
1000 times, in order to ensure that the full range in con- 
tinuum variability behaviour is sampled and to allow the 
construction of probability distribution functions in each of 
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Figure 13. As for figure fT2l for Mg II (hashed area) and C IV. 



the desired quantities that we wish to measure. For each 
continuum - emission-line light-curve combination we have 
calculated the peak and centroid of the cross-correlation 
function ffigures [T2lll3j) . Additionally, from the velocity re- 
solved time-variable light-curve, we have constructed mean 
and root-mean square profiles from which measurements of 
their full width at half maximum (fwhm) and line disper- 
sion (ai) have been determined. Thus, from each simulation 
we have two estimates of the continuum-emission-line delay 
(or characteristic size of the BLR) , and four estimates of the 
velocity dispersion (two from the mean profile and two from 
the rms profile). From these we construct eight estimates of 
the BLR mass in the standard fashion using equation [6] 

For a few individual objects, intensive ground and 
space-based spectroscopic monitoring campaigns have been 
employed to determine black holes masses under the as- 
sumption of virialised gas motion for broad emission-lines 
spanning a broad range in ionisation state (e.g. Clavel 
et al. 1991, Ulrich et al. 1991, Peterson et al. 1991, Krolik 
et al. 1991). For NGC 5548 Peterson and Wandel (1999) 
showed that while lines of differing ionisation state dis- 
play differences in both their line dispersion and response 
timescale they all appear to follow (more or less) the same 
virial relation. Furthermore, Peterson et al. 2004, and Collin 
et al. 2006, using 13 years of monitoring data of the op- 
tical emission-lines in NGC 5548 showed that differences 
in the measured lags and line widths from one season to 
the next yielded virial products which are effectively con- 
stant within the errors, demonstrating that the reverbera- 



tion masses while not necessarily accurate, are at least ro- 
bust (reproducible). However, as noted by these authors, it 
is not immediately obvious which of the measured quantities 
(CCF peak or centroid, fwhm, ai of the mean or rms pro- 
file) yield the most accurate black hole masses, nor which of 
the emission-lines should be used in mass determinations. 

Theoretical considerations suggest that measurements 
from rms profiles rather than mean profiles should yield 
more accurate estimates of the virial product, since they iso- 
late the variable part of the emission-line from non-variable 
components (Fromerth and Melia 2000). However, because 
of the lower S/N of rms profiles, a consequence of the small 
variability amplitude, measurements from rms profiles tend 
to be more uncertain. Collin et al. (2006) argued that pro- 
vided the non- varying components can be isolated from the 
mean profile, then the mean profile shoul d y ield estimates 
of the virial product of a similar precision 16 !. Additionally, 
Peterson et al. (2004) showed that the line dispersion ai can 
be measured with greater precision than the line fwhm and 
provides a better match to the virial relation, which sug- 
gests that estimates of the virial product based on ai should 
have smaller uncertainty. Finally, the centroid of the CCF, 
equivalent to the centroid of the 1-d response function (or 
the emissivity-weighted radius for a linear response) is the 
generally preferred quantity used for estimating the 'size' of 
the BLR since it is less sensitive to the inner BLR radius 
than is the peak of the CCF (Peterson et al. 1998), but see 
Welsh 1999, for an alternative viewpoint. Indeed, for ob- 
jects in which monitoring campaigns have been performed 
for multiple lines, Peterson et al. 2004 found that the tight- 
est virial relation was found for virial products determined 
from measurements of the centroid of the CCF (r cen t) and 
the line dispersion (ai) of the rms spectrum. 

In summary, from an observational perspective virial 
products based on measurement of the CCF centroid and the 
dispersion of the rms profile remain the quantities of choice. 
However, it seems at least plausible that the choice of these 
quantities may simply reflect limitations in the observational 
data, be it signal-noise or sampling. Here we attempt to 
address the question as to which of the measured quantities 
and for which lines yield the most accurate black hole masses 
from a modelling perspective. 



5.3.1 CCF centroid or peak 

In terms of the measured delay, or characteristic size of the 
BLR, it has long been known that the CCF peak (or lag) is 
a less reliable measure of the characteristic BLR size than 
is the centroid, and tends to be biased towards the inner 
BLR radius (see e.g. Edelson and Krolik 1988; Perez, Robin- 
son and de la Fuente 1992; Melnikov and Shevchenko 2008) 
though see Welsh 1999 for an alternative point of view. By 
contrast, the centroid of the CCF is directly related to the 
centroid of the response function (the responsivity weighted 
radius, Koratkar and Gaskell 1991) and thus should be a 
more accurate representation of the 'size' of the region re- 
sponding to continuum variations (the correspondence is not 



16 Collin et al. 2006, as with previous authors find that mean 
profiles are typically 20% broader than rms profiles, possibly due 
to the presence of low responsivity gas in the inner BLR. 
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exact even for the centroid and depends strongly on the con- 
tinuum variability behaviour exhibited during the observing 
campaign, see e.g. Perez, Robinson and de la Fuente 1992). 

Figures [T2l413l show the distribution in the measured de- 
lay in terms of both the centroid and peak delay (or lag) for 
each line over a range of BLR inclinations. Also shown is the 
centroid of the 1-d responsivity-weighted response function 
(vertical dotted line for He n and C IV, vertical dashed line 
for H/3 and Mg n) for the steady-state for the same inclina- 
tion. For low inclinations the CCF centroid for all lines bar 
He II tracks the centroid of the response function reason- 
ably well, while for larger inclinations (i > 10 degrees), the 
measured centroid displays larger scatter and is systemati- 
cally smaller than the true centroid in all lines. The discrep- 
ancy between the centroid of the 1-d responsivity-weighted 
response function and the centroid of the CCF is notably 
worse at large inclinations for He n and C IV, which in our 
model preferentially form at small BLR radii. Furthermore 
while at low inclinations the peak of the CCF also tracks the 
centroid of the response function reasonably well, it performs 
far worse at larger inclinations. For He n the peak delay is 
systematically smaller than the centroid of the 1-d response 
function regardless of inclination. Our adopted continuum 
variability timescale of 40 days is close to the centroid of the 
1-d response function for H/3 and Mg n respectively and is 
the most likely explanation for their improved performance 
in tracking the continuum variations when compared to C IV 
and He II, particularly at large inclinations. 
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5.3.2 Mean or rms profiles 

Figures [T4HT51 show the derived black hole mass distribution 
functions (8 independent estimates) as a function of line of 
sight inclination i for all four lines. Because the mass dis- 
tribution functions are typically skew the points and their 
errors here represent the median and 68% confidence lim- 
its respectively. The input black hole mass is indicated by 
the horizontal dashed line. As expected, the derived Mbh 
for all four lines shows a strong dependence on inclination, 
with black hole mass estimates systematically lower than the 
input mass at low line of sight inclinations. Mass estimates 
based on measurements using the line fwhm (from both the 
mean and rms profile) show a greater dependence on incli- 
nation at small inclinations, than do mass estimates based 
on measurements of the line dispersion ai (ie. as noted by 
Collin et al. 2006, ai is a less biased indicator of mass). 

For our model, fwhm/ai « 1 at i=2 degrees, and in- 
creases to a maximum for inclinations i > 20 degrees after 
which time fwhm/ai is approximately constant (see also 
figure [TO]) . Thus the line shape is largely insensitive to incli- 
nation effects at large inclinations but is a strong function 
of inclination at small inclinations (i < 20 degrees), inclina- 
tions for which scale-height dependent turbulence may sig- 
nificantly modify the line profile (see § 15. 1|) . The origin of 
this effect in our model is due entirely to the implementa- 
tion of Gravitational redshift and transverse Doppler shift 
(see ^5.1|) . In the absence of these effects, only the line width 
(and not line shape) is a strong function of inclination. 

Since fwhm/ai > 1 for all lines (except at the lowest 
inclinations) mass estimates based on measurements of the 
fwhm are generally larger than those based on ai . Black hole 
mass estimates based on the line fwhm systematically un- 



Figure 14. 8 estimates of the black hole mass, Mbh (median, 
and lcr confidence interval) based on the virial product Rbi r 5v 2 /G 
as a function of inclination i for our fiducial BLR model (see Table 
1), and based on 1000 simulated light-curves at each inclination. 
Our calculated virial products (Eqtn. 6) employ two measures of 
Sv, fwhm and ai, and two measures of Rbw (= CT )-> the CCF 
centroid and CCF peak. Together these provide 4 estimates of 
Mbh for both the mean and rms profiles (giving 8 estimates in 
total). Individual panels show results for H/3 (open circles, solid 
black line), and He II A4686 (open squares, red dashed line). 



derestimate the mass at low line of sight inclinations, while 
systematically overestimating the mass for line of sight in- 
clination i > 20 degrees (see figures [LllJ15[) . Mass estimates 
based on measurements of the line dispersion, ai, systemat- 
ically under-predict the black hole mass at all inclinations 
and for all lines, with the largest discrepancies at the low- 
est inclinations. The larger scatter in measurements of the 
fwhm and ai of the rms profile is as expected (a reverber- 
ation effect), but does not necessarily lead to a larger scat- 
ter in the derived mass, if variations in the line widths are 
compensated for by corresponding changes in the measured 
delay. 

The use of CCF centroid or peak in the virial relation 
in general makes only a small difference to the estimated 
masses (see e.g. left-hand or right-hand panels of figures [12]- 
I13|) though the exact behaviour depends both on the line in 
question and details of the continuum variability behaviour 
during a single model run. For example, for He n the CCF 
centroid is on average larger than the CCF peak at all in- 
clinations, while for C IV, the CCF centroid is on average 
larger than the CCF peak only for i > 10 degrees. Similarly, 
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Figure 15. As for figure [14] but now showing the results for 
Mg II A2798 (open circles, solid black line), and C IV A1550 (open 
squares, red dashed line). 

for H/3, the CCF centroid is on average larger than the CCF 
peak for inclinations in the range 10 < i < 50 degrees, and 
smaller than the CCF peak otherwise. For Mg n the CCF 
centroid is on average larger than the CCF peak for inclina- 
tions in the range 5 < i < 30 degrees, and smaller than the 
CCF peak otherwise. Thus when considering all lines, we 
expect to find on average larger black hole masses at large 
inclinations when using the CCF centroid, while at low incli- 
nations the CCF centroid yields lower mass estimates than 
the CCF peak for all lines except He II. 

We have recast the derived Mbh from our reverbera- 
tion simulations of our fiducial BLR model to illustrate the 
variation in virial scale factor / as a function of inclination 
for each line. Here we restrict ourselves to four estimates 
only, the four used by Collin et al. 2006, in which velocity 
dispersion measurements (fwhm, o~i) have been taken from 
the mean and rms profiles, and the centroid of the CCF 
is used as a proxy for BLR size. The results of this exer- 
cise are displayed in figure [TBI where the points and their 
errors once again indicate the median and 68% confidence 
limits respectively. Increasing the characteristic variability 
timescale for the continuum variations from 40 to 80 days 
only marginally alters the virial scale factor for H/3 (smaller 
for the line dispersion and larger for the line fwhm). 

Continuum variability blurs the differences between the 
predicted virial scale factors for each line to the point at 
which they are virtually indistinguishable within the errors 
at inclinations larger than ^20 degrees. Thus when measure- 



Figure 16. Four estimates of the virial scaling factor f, for four 
different lines driven by 1000 continuum light-curves with a char- 
acteristic variability time scale r c h a r °f 40 days. Black filled cir- 
cles (H/3), blue triangles (Mg il), green squares (C iv), and cyan 
crosses (He il). The horizontal solid and dashed lines represent 
the average H/3 scale factors, /(H/3) and their uncertainties, as 
derived for the Onken et al. sample by Collin et al. 2006. 

ment errors are included we predict little if any difference be- 
tween the f -factors (or mass estimates) derived for different 
lines. Importantly, /-factors derived from the fwhm of the 
rms and mean profile show significantly less scatter and a 
closer correspondence between the lines than those based on 
ai, particularly at low inclinations. This consistency between 
mass estimates for different emission lines derived using the 
line fwhm has been confirmed observationally (Assef et al. 
2011). Because our light curves are well-sampled, there is 
only a small spread in the measured values of fwhm and 
gi of the mean profile. Thus the dispersion in the computed 
virial scale factors from measurements of mean profiles is 
mostly due to variations in the location of the centroid of 
the CCF. For light-curves which are only sparsely sampled, 
we would expect larger variations in fwhm and o~i in the 
mean profiles, though still smaller than those measured for 
rms profiles. Figure [16] shows that for mass estimates based 
on measurements of the line fwhm, the /-factors for all 
lines are consistent within the errors (/ db Sf) with those 
reported by Collin et al. (2006) for inclinations in the range 
20 < i < 35 degrees. For estimates based on measurements 
of the line dispersion ai, the acceptable range narrows to 
20 < i < 30. 

The average virial scale factor for the Decarli et al. sam- 
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pie is f d = 1.6 ± 1.1 for H/5 and / d = 2.4 ± 0.16 for C iv, 
that is /-factors tend to be larger for C IV than H/3. For 
our fiducial model, in the presence of turbulence f(H/3)(ai) 
< f(C iv)(ai) for i < 20 degrees for both the mean and 
rms profile (figure [16]) , because of the stronger dependence 
of the H/3 line width on turbulence due to it forming at pref- 
erentially greater distances and scale height. By contrast f- 
factors derived from the fwhm of the mean and rms profile 
indicate f(C w)(fwhm)~ f(Hf3)(fwhm), though the cor- 
respondence is generally worse for the lowest inclinations. 

The larger scatter in the /-factor between the low and 
high ionisation lines at low i derived using a\ (figure [16] - 
left hand panels) highlights the reduction in sensitivity of 
the line-dispersion to inclination in the presence of strong 
turbulence (which mainly effects the low ionisation lines). 
We suggest that objects for which the lines show large dif- 
ferences in their virial scale factors (as derived using the line 
dispersion) may be indicative of low inclination sources. 



6 DISCUSSION 

A bowl-shaped BLR geometry provides an elegant solution 
to the smaller than predicted dust reverberation sizes by 
decreasing the measured delays without altering the dust 
formation radius. At the same time, since material now lies 
away from the observers line of sight, they readily reproduce 
the observed lack of response on short timescales evident 
in the 1-d and 2-d response functions of the strong optical 
recombination lines. 

Differences in the radial surface emissivity distribution 
and line anisotropy among the low and high ionisation lines 
in the context of a bowl-shaped BLR, results in large dif- 
ferences in the form of the 1-d response functions at fixed 
inclination (see figure A3, appendix). For geometrically thin 
discs, the characteristic variability timescale as determined 
from the centroid of their 1-d response functions is inde- 
pendent of inclination for isotropically emitting clouds, and 
increases with inclination otherwise. For bowl-shaped ge- 
ometries, the centroid of the 1-d response function increases 
with inclination even for isotropic emission. Moreover for 
the radiation pattern adopted here, at small inclinations, 
increased line anisotropy can reduce the centroid of the re- 
sponse function relative to the isotropic case due to the re- 
duced contribution of gas at large elevations which lies closer 
to the line-of-sight. 

While the LILs show an absence of response on short 
timescales in their 1-d response function, the 1-d response 
functions for the HILs, which are formed at smaller BLR in 
a more flattened distribution, resemble those of disc-shaped 
BLR geometries, with a significant response even on short 
timescales (though the response does decline to zero at zero 
delay). Thus this model reproduces the observed differences 
in the location of the peak response among low and high ioni- 
sation lines reported in the literature (e.g., Krolik et al. 1991, 
their figure 10 and 11). We note that here we have assumed 
a locally linear line response and that the precise form of 
the 1-d response function may be modified in the event of 
significant non-linear effects (this may include a luminosity- 
dependent continuum shape as well as the incident contin- 
uum flux-dependent effects already mentioned). 

The additional effects of GR and TDS enhance the red- 



wing response producing line profiles with extended red- 
wings, as are sometimes observed in type 1 AGN (Kol- 
latschny 2003) . The strength of the red-blue asymmetry de- 
pends primarily on the line formation radius, being stronger 
for lines formed at small BLR radii (i.e., steep radial emis- 
sivity distributions). GR and TDS together provide signif- 
icant line width (~ several hundred km s _1 ) even for face- 
on geometries with pure-planar Keplerian motion. More im- 
portantly, GR and TDS introduce a strong inclination de- 
pendence to the line profile shape at low inclinations. In 
the absence of these effects the shape of the emission-line 
profiles are independent of inclination for flattened/bowl- 
shaped BLR geometries assuming isotropic emission, though 
of course their line widths will show a strong inclination de- 
pendence (the observed line of sight velocity, v bs, varies as 
vsini). 

Mass determinations for flattened/bowl-shaped BLR 
geometries show a strong dependence on inclination particu- 
larly at small inclinations where the mass estimates are sys- 
tematically smaller than the input model. In general, when 
reverberation effects are included, mass estimates based on 
measurements of the fwhm underestimate the mass at low 
inclinations, and overestimate the mass at high inclinations, 
while those based on measurements of the line dispersion o~i , 
systematically underestimate the mass at all inclinations. 
Mass estimates based on measurements of the emission line 
fwhm (rms or mean profile) are larger since in general 
fwhm/ai > 1. Our model also confirms the result of Collin 
et al. 2006, showing that mass estimates based on measure- 
ments of gi are less biased than those determined from mea- 
surements of the fwhm, because of the weaker dependence 
of a i on inclination, particularly at low inclinations. For our 
simulations, we find a better correspondence between mass 
determinations derived from different lines particularly at 
low inclinations, if measurements of the virial product are 
performed using the fwhm of the mean or rms profile. The 
form of emission line anisotropy adopted in this model leaves 
the profile unchanged, but causes the characteristic response 
timescale for the line to increase with increasing inclination. 
This effect is smaller than the corresponding increase in ve- 
locity with inclination and consequently both Mbh and / 
are only weakly dependent on line anisotropy (e.g. figure [TT1 
lower left panel). 

Turbulence, as implemented here, modifies the shape of 
the 2-d response function and emission-line profile, by mov- 
ing lower velocity gas that responds on long timescales to 
larger line of sight velocities (hence broadening the lines). 
The overall effect on line shape is line dependent, and is 
largely determined by the line formation radius, so that 
while turbulence may increase the ratio fwhm/ai for lines 
formed at small BLR radii (small scale heights) , the general 
effect is to reduce fwhm/ai for lines formed at large BLR 
radii (large scale heights), so that at low line of sight incli- 
nations, the line profiles are characterised by narrow cores 
and extended line wings (i.e., Lorentzian), similar to those 
seen in NLSyls, see figure 8 and 35.11 for a full discussion 
of these effects. Because turbulence randomises the velocity 
field, turbulence acts to reduce the i>sinz dependence of the 
line dispersion in flattened BLR geometries. 

To summarise, in the absence of turbulence, emission- 
lines with steeper emissivity distributions yield : (i) larger 
estimates for the central black hole masses, (ii) smaller virial 
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Figure 17. fwhm/ai(Hf3) versus (i) cr^(H/3) (top panels), and 
(ii) fwhm(YLfi) (lower panels) for the rms (left hand panels) and 
mean profiles (right hand panels) for our fiducial bowl-shaped 
model with turbulence parameter 6 tur b = 2, driven by 1000 sim- 
ulated continuum light-curves, and observed at line of sight in- 
clinations 2 (yellow open circles), 5 (green), 10 (cyan), 20 (blue), 
30 (magenta), 40 (red), and 50 (black) degrees. At a given in- 
clination the rms profile shows a larger range in fwhm/ai, a 
consequence of reverberation within the spatially extended BLR. 
A far smaller range is seen for fwhm/ai at fixed inclination for 
the mean profile, as measured from the same simulations (same 
number of points). For i > 20 degrees, the profile shape is largely 
independent of inclination. 



scale factors (f- values) and (iii) smaller fwhm/a (see e.g. 
figure [TO]) . 



6.1 Line shape as an orientation indicator 

Collin et al. 2006, showed that the shape (fwhm/ai) of 
the broad H/3 line in reverberation mapped AGN varies 
by a factor of a few, and tends to be smaller in Narrow 
Line objects (fwhm/ai < 2.35) than in Broad Line objects 
(fwhm/ai > 2.35), where fwhm/ai = 2.35 is the value ap- 
propriate for a Gaussian profile. The boundary separating 
the narrow line and broad line objects, hereafter population 
1 and population 2 sources, roughly corresponds to a ai of 
2000 km/s, and is broadly similar to the division of AGN 
into population A population B sources by Sulentic et al. 



(20000- In the latter scheme, population A sources have 
profiles characterised by narrow cores and broad wings (ie. 
are more Lorentzian in shape) while population B sources 
appear more flat-topped (i.e., more Gaussian). 

By estimating the virial product using fwhm and ai of 
the mean and rms profiles of each sample and comparing to 
the Mbh values obtained using stellar velocity dispersion 
measurements for the same AGN, Collin et al. (2006) cal- 
culated virial scale factors for both Pop 1/A and Pop 2/B 
samples, as well as the sample as a whole. A key result of 
their work is that the virial scale factor determined from 
measurements of the line dispersion (mean or rms profile) is 
independent of line shape, or AGN type. By contrast, virial 
scale factors based on the line fwhm show significant differ- 
ences among the different AGN populations (pop 1/A, pop 
2/B) and is thus a more biased estimator of the velocity 
dispersion. 

Collin et al. (2006) speculate that the factor of w 3 dif- 
ference in the virial scale factor derived using the fwhm for 
population 1/A and population 2/B sources may be due to 
the increased sensitivity of the fwhm to inclination effects 
or Eddington ratio, since population 1/A sources tend to be 
high Eddington rate sources. They suggest that the rapid 
decrease in fwhm/ai at small ai arises in a 2-component 
BLR, comprising a disc component producing the emission- 
line core, and an isotropic (possibly a wind) component pro- 
ducing the line wings. Since the core of the emission-line is 
more sensitive to the fwhm and arises in a flattened con- 
figuration, it will show a strong dependence on inclination. 
Conversely, the wings of the line are more sensitive to ai, 
and arise in an isotropic component which is inclination in- 
dependent. We here offer an alternative explanation. The 
rapid decrease in fwhm/ai at small o\ arises naturally in 
planar disc-like (thin or thick) geometries in which the ef- 
fects of GR and TDS are taken into consideration. At small 
inclinations ai is independent of inclination while the ratio 
fwhm/ai increases rapidly (factor of ~ 5). That is, the line 
fwhm is more sensitive to changes in inclination for small 
inclination angles than is the line dispersion. This effect can 
be seen more easily in figure [17] where we show the ratio 
fwhm/ai for H/3, as measured from the mean and rms pro- 
files (left and right-hand panels respectively) as a function of 
gi (upper panels) and fwhm (lower panels). At low inclina- 
tions i < 20 degrees, ai is effectively constant. Variations in 
line shape are then entirely due to changes in the line fwhm 
(compare the upper and lower left panels of figure [T7)) . At 
higher inclinations, the line shape is nearly independent of 
inclination. This trend is seen in measurements from both 
the mean and rms profiles. The reduced scatter in the ra- 
tio fwhm/ai for the mean profile relative to the rms profile, 
suggests that measurements taken from the mean profile are 
less sensitive to reverberation effects with the spatially ex- 
tended BLR. Turbulence shifts all velocities to higher values 
while its effect on fwhm/ai depends on the scale height at 
which the line forms. For large turbulent values fwhm/ai 



17 Decarli et al. (2008) claim that the correlation found between 
fwhm/ai and a\ is an artifact of the fitting procedure employed 
to quantify the velocity dispersion, though as we have already 
pointed out differences in fwhm/ai are to be expected in BLR 
geometries bounded by an inner and outer radius. 



22 Goad, Korista and Ruff 



tends toward 2.35 (a Gaussian profile), as appropriate for a 
randomised velocity field. In any event, it is difficult to see 
how the bottom right hand corner of figure [TTl can be pop- 
ulated in this model, in agreement with observations (see 
Collin et al. 2006, their figure 3; Kollatschny & Zetzl 2011, 
figures 1-3; Peterson 2011). Remarkably, for our model the 
break toward smaller fwhm/ai occurs at a <ji of between 
1500 and 2000 km/s, similar to the boundary between pop- 
ulation 1 and population 2 sources defined by Collin et al. 
2006 (see also Peterson 2011). Figure 17 shows the range 
in profile shape as a function of fwhm and ai for a single 
source with fixed Mbh, viewed at a range of line of sight 
inclinations and in the presence of reverberation effects (ie. 
approximating the results of many observing seasons). By 
comparison, Figure 3 of Collin et al. (2006) indicates the 
mean profile shape observed in AGN which differ in both 
Mbh and line of sight inclination. Yet our figure 17 and fig- 
ure 3 of Collin et al. are remarkably similar in appearance. 
Increased mass will tend to shift points in the upper left 
panel of Figure 17 toward higher velocity dispersion, effec- 
tively filling in the left hand side, while leaving the line pro- 
file shape at larger velocity dispersions, virtually unchanged, 
as is observed. Similarly, the lower left panel of Figure 17 
bares a striking resemblance to Figure 1 of Kollatschny and 
Zetzl (2011). Yet in their model a Lorentzian component 
was input by hand, while here it arises naturally in low in- 
clination systems for lines formed at large BLR radii and in 
the presence of scale height dependent turbulence. 

As expected reverberation effects can alter the emis- 
sion line profile shape, though for our simulations the range 
in fwhm/ai at a fixed inclination is by comparison with 
NGC 5548, relatively modest (-10-20%). We suggest that 
the far larger range in fwhm/ai observed over the 13 
year monitoring campaign of NGC 5548, which at times 
crosses the boundary between population 1 and population 
2 sources, is a consequence of changes in the radial surface 
line emissivity distribution in response to changes in the ion- 
ising continuum flux. We also find that at a fixed inclination 
the range in fwhm as determined from the rms spectrum is 
far larger than that for ai (figure [TTl left panels) , suggesting 
that in our model, the line dispersion is also less sensitive to 
reverberation effects within the physically extended BLR. 
The far smaller spread in fwhm and ai for the mean spec- 
tra (figure \T7\ - right panels) is mainly due to the absence 
of contaminating non- variable components and the fact that 
the light-curve is uniformly sampled with no gaps between 
observations. 

6.2 Anomalous narrow-line quasars 

In a recent study on anomalous narrow- line quasars from 
the Sloan digital sky survey, Steinhardt & Silverman 2011 
(arXiv:1109.1554 ) showed that the virial masses determined 
using the familiar scaling relations derived for the broad H/3 
and Mg n emission- lines differ by up to 0.5 dex in these ob- 
jects. ANLs are identified by a broadened narrow H/3 line 
relative to Mg II. This broadened narrow- line component is 
found to be well- correlated with broad H/3 in these objects. 
Since H/3 and Mg n are expected to form under similar phys- 
ical conditions, one would normally expect a similar relation 
to be evident in Mg II. Steinhardt & Silverman 2011 argue 
that the absence of such a correlation in Mg II may cast 



doubt as to the validity of virial mass estimates based on 
measurements of the H/3 line width, arguing that some of 
the broadening may be due to the presence of a wind. Ther- 
mal emission from a wind lacking an ionisation front would 
produce an additional H/3 component without contributing 
significantly to the Mg n emission. 

Bowl-shaped geometries may provide an alternative so- 
lution. If as expected Mg n forms at large radial distances, 
and turbulence is significant at large scale heights, then one 
can envisage a situation in which the Mg n emission is dom- 
inated by a turbulent component. The Mg n line width 
would then be largely independent of inclination effects. By 
contrast if H/3 emission arises at smaller radial distances 
(smaller scale heights), the lower turbulence and more flat- 
tened spatial distribution would ensure a strong inclination 
dependence even at modest inclinations (GR and TDS are 
less important for these lines due to their larger formation 
radii). Thus, ANLs may represent near face-on objects with 
a significant turbulent component in the Mg n forming re- 
gion. 



7 SUMMARY 

We have explored the observational characteristics of a class 
of broad emission line region geometries in which the line 
emitting gas spans the region between the outer accretion 
disc and the inner edge of the dusty torus, by occupying 
an effective (though not necessarily smooth or continuous) 
surface the scale height of which increases with increasing 
radial distance (ie. a bowl-shaped BLR geometry), similar to 
the configuration first proposed by Gaskell (2009). A Type 1 
AGN spectrum is then observed for line-of-sight viewing an- 
gles which peer over the rim of the bowl. Such a configu- 
ration provides quite naturally the necessary high covering 
fractions for the observed emission-line strengths, without 
obscuring the view of the central engine (for a reasonable 
range in line-of-sight viewing angles). The gas dynamics are 
here assumed to be dominated by gravity, and we include 
in our model the effects of gravitational redshift, transverse 
Doppler shift and scale-height dependent turbulence. While 
we do not exclude possible contributions from infall, out- 
flow, or line-driven wind contributions to the BLR kinemat- 
ics within individual AGN, we consider here what is most 
likely the underlying and dominant effect of a gravitationally 
bound BLR velocity-field consisting of circularised orbital 
motion, as might be expected if there are significant dissipa- 
tive forces acting within the BLR. Importantly, we introduce 
a macroscopic turbulent component to the cloud motion 
which imparts substantial scale height to the BLR at large 
radii, allowing it to intercept sufficient ionising continuum 
radiation to explain the observed emission-line strengths. 

By breaking spherical symmetry, a bowl-shaped geome- 
try simultaneously provides a simple solution to the shorter 
than predicted dust reverberation timescale for the inner 
dusty torus and the absence of a significant response of the 
broad optical recombination lines on short timescales. Grav- 
itational redshift and transverse Doppler shift offer a mecha- 
nism for providing substantial line width, even in the case of 
purely transverse motion. They also introduce a strong red- 
blue asymmetry into the 2-d response function and emission- 
line profile in the form of an enhanced red-wing response at 
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short time-delays similar to that seen in the recently recov- 
ered 2-d emission-line response function and emission-line 
profile of the NLSyl Arp 151 (Bentz et al. 2010, Brewer 
et al. 2011). These effects are most important for lines with 
steeper radial emissivity distribution which therefore form at 
small BLR radii (the HILs, e.g. He II, C iv) and for those sys- 
tems viewed at low line of sight inclination. Differences in the 
line formation radius between the high- and low-ionisation 
lines lead to strong differences in the shapes of their 2-d 
and 1-d emission-line response functions and emission-line 
profiles. 

While the LILs show a deficit of response on short- 
timescales, the HILs, which form toward the base of the bowl 
in a more flattened configuration, show significant response 
on short-timescales, and their 2-d and 1-d response func- 
tions appear similar to those found for geometrically thin 
discs. Additionally, GR and TDS for flattened BLR geome- 
tries (ie. bowl-shaped, thick or thin disc geometries), result 
in line profile shapes which display a strong inclination de- 
pendence at low line of sight inclinations, with fwhm/ai 
decreasing with decreasing i. At larger inclinations, the line 
profile shape is nearly independent of inclination. Thus in- 
clination, in the presence of GR and TDS may provide a 
natural explanation for the observed differences in line pro- 
file shapes among AGN in general and importantly may also 
explain the observed differences in the line profile shapes 
used to discriminate between pop A and pop B sources (Su- 
lentic et al. 2000, Collin et al. 2006). We suggest that some 
of the systems previously identified as high Eddington rate 
sources based on their profile shape, may simply be type 
1 objects viewed at low line of sight inclinations. We note 
that the break between inclination dependent and inclina- 
tion independent line profile shapes (ai ~ 2000 km/s, see 
figure fT7|) is remarkably similar to the boundary separating 
the high- (pop A) and low- (pop B) Eddington rate sources. 
Also, we speculate that high Eddington rate sources may 
have TORs with larger than typical covering fractions, ne- 
cessitating smaller observed viewing angles i for un-obscured 
lines of sight. 

For a bowl-shaped BLR with significant scale-height de- 
pendent turbulence, we expect the low-ionisation lines (e.g. 
H/3, Mg n) to display line shapes which are more Gaussian in 
shape (larger fwhm/ai) than the high ionisation lines (e.g. 
C IV, He n), as is observed (Decarli et al. 2008). This arises 
because the low-ionisation lines form at larger BLR radii 
(and hence larger scale heights) where the turbulent (ran- 
dom) contribution to the velocity field is larger. By contrast 
the high-ionisation lines form at small BLR radii (near the 
base of the bowl) where the turbulent component is sub- 
stantially reduced (for example, compare the fwhm/ai as 
a function of inclination i for the low- and high-ionisation 
lines with and without turbulence, figure [TOl lower panels). 
As can be seen, as well as significantly broadening the line 
profile, turbulence softens the strong dependence of line pro- 
file shape on inclination at low inclinations for the low- 
ionisation lines. By moving low velocity gas (at large BLR 
radii) to larger line of sight velocities, turbulence removes 
the strong horns and shoulders characteristic of pure planar 
Keplerian motion while reducing any delays that may be 
present between the wings and core of the line. For low line 
of sight inclinations, the turbulent contribution to the line 
of sight velocity is larger than the Keplerian contribution 



(for motion confined to the disc mid-plane) and turbulent 
broadening acts to produce emission line profiles which are 
Lorentzian in form (fwhm/ai ~ 1, e.g. figure [8]), with ex- 
tended broad wings and narrow cores. As far as we aware, 
this is the first study to provide physical motivation for the 
fitting of emission-line profiles using Lorentzians. A signifi- 
cant scale-height dependent turbulent component may help 
explain the reported weakness in the correlation between 
the measured widths of the low- and high- ionisation lines 
(e.g., H/3 and C iv). In cases of sufficiently large contribu- 
tions from this component, H/3 might even be broader than 
C IV, a finding which when taken at face value, appears at 
odds with what one would normally expect based on pho- 
toionisation calculations of a gravitationally bound BLR. 

A strong motivation for this work was to try to iden- 
tify which lines and which of the measured quantities yield 
the most accurate and most reproducible mass estimates for 
Mbh- The time-stationary model indicates that : 

(i) Different emission lines predict different virial mass 
estimates and hence different virial scale factors (figure [TO)) . 

(ii) Mass estimates made from measurement of the emis- 
sion line fwhm show a strong inclination dependence, lead- 
ing to an underestimate of Mbh at low inclinations, and an 
overestimate of Mbh at large inclination. The largest dis- 
crepancies are for low line of sight inclinations and in the 
absence of turbulence for lines formed at large BLR radii 
(ie. the LILs, see figure fTQ|) . 

(iii) Mass estimates made from measurement of the emis- 
sion line dispersion ai, underestimate Mbh at all inclina- 
tions, but show a much weaker dependence on i than the 
fwhm (see also figures U4H15p . confirming the results of 
Collin et al. (2006), that mass estimates based on ai are 
less biased. 

(iv) Measurements of the emission line fwhm (mean or 
rms profile) predict larger Mbh than measurements based 
on ai, since in general fwhm/ai > 1, except at the smallest 
inclinations. 

(v) In the presence of turbulence, black hole mass esti- 
mates based on measurements of emission lines formed at 
large BLR radii (ie. H/3, Mg n), are less sensitive to inclina- 
tion effects than those formed at small BLR radii (ie. C IV, 
He n), and for which the turbulent contribution to the ve- 
locity field is weaker. 

(vi) In the absence of turbulence, emission lines formed at 
small BLR radii (ie. the HILs) yield larger Mbh (figure [TO]) 
and therefore smaller virial scale factors. 

When the effects of emission-line variability are taken into 
account, we find that: 

(i) The cross correlation function (CCF) centroid tracks 
the continuum variations more closely than the CCF peak, 
and shows a better correspondence with the centroid of the 
emission line response function in lines whose variability 
timescale is more closely matched to the characteristic vari- 
ability timescale of the driving continuum (ie. the LILs, fig- 
ures [12HT3} . However, the choice of CCF centroid over CCF 
peak has little overall effect on the predicted mass. 

(ii) The best correspondence between Mbh estimates 
(and hence /-factors) for different lines, arises from mea- 
surements of the fwhm of the mean or rms profile, consis- 
tent with results reported in the literature from observations 
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(e.g. Assef et al. 2011). The stability of the mean profile sug- 
gests that measurement of the fwhm of the mean profile is 
the preferred estimate of Mbh (figure [16]). 

(iii) The largest discrepancy between the derived virial 
scale factors between different emission lines arise from mea- 
surements of <ji in low inclination systems (i < 20 de- 
grees, figure 16). We suggest that some of the systems pre- 
viously identified as high Eddington rate sources based on 
the shapes of their low-ionisation emission lines, may simply 
be more typical type 1 objects viewed at low line of sight 
inclinations. 

Our findings at first glance appear at odds with those 
presented in Peterson et al. (2004) who found that the tight- 
est relation between the virial products derived for different 
lines in a single source is found for those estimated using the 
CCF centroid, and the line dispersion of the rms spectrum. 
Note that the difference here is that we calculate the time- 
delay (centroid or peak) of the CCF and line width (ai or 
fwhm) of the mean and rms profile for 1000 realisations of 
the driving continuum light-curve, and not the variation in 
centroid and line width within a single realisation. Due to 
the long duration of our light-curves, any biases introduced 
by windowing effects are minimised. In addition, the model 
presented here utilises a locally linear response approxima- 
tion. Gross changes in the radial surface line emissivity dis- 
tribution and in the inner and outer radius in response to 
large variations in the ionising continuum flux variations will 
be treated elsewhere (Goad and Korista, in prep.). 

A comparison between the H/3 virial scale factors pre- 
sented in our simulations and the empirically estimated 
virial scale factors based on the AGN with reverberation 
data presented in Collin et al. (2006) suggests a rather nar- 
row range of typical inclination angles for Type 1 objects, 
i « 20-35 degrees or so, in agreement with the expected 
range found via other techniques. 
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APPENDIX A: RESPONSE FUNCTIONS AND 
EMISSION-LINE PROFILES FROM OUR 
SIMULATIONS 

In figures IA1HA3I we show the model 2-d and 1- 
d responsivity-weighted response functions and variable 
emission-line profiles for our fiducial BLR geometry for each 
of the four lines described in the text and line of sight view- 
ing angles in the range 2-40 degrees. The models include the 
effects of TDS and GR and have been calculated assuming 
a turbulence parameter 6 t urb = 2. 



There are notable differences between the 2-d and 

1- d responsivity-weighted response functions and variable 
emission-line profiles for all of the lines. Broadly speaking 
the lines fall into two categories with the dividing line set 
by the scale height at which the line forms. H/3 and Mg n 
are formed at large scale heights. Consequently at low incli- 
nations they display broader 2-d response functions at large 
time-delays and Lorentzian profiles. The effects of TDS and 
GR are less prominent for these lines and the line profiles 
appear more symmetric. The 2-d response functions at low 
inclination indicate a large deficit in response at short time- 
delays. The peak in the response function is located at the 
outer edge on the side of the bowl nearest the observer, and 
moves to smaller delays as the inclination increases. At small 
inclinations, the turbulence is large enough to fill in the gap 
between the horns normally present in the emission-line pro- 
files for bowl-shaped geometries. At larger inclinations, the 
horns re-appear because of their large separation in velocity 
space. Cranking up the turbulence would once again remove 
the peaks from the emission-line profile. 

C IV and He n have steeper emissivity distributions and 
consequently form at smaller BLR radii, which in our model 
implies small scale heights. Since these lines originate in a 
more flattened distribution, the effects of turbulence on the 

2- d and 1-d response functions are less significant. The 2-d 
response functions therefore display a strong red-blue asym- 
metry on short time-delays with an enhanced red-wing re- 
sponse. Their emission-line profiles are broader in compari- 
son to H/3 and Mg n at large inclinations, and the enhanced 
red-wing response is still evident even at large inclinations. 
Large differences are also notable in the 1-d response func- 
tions, where at large inclinations the response functions of 
C IV and He n show a strong resemblance to those found 
for geometrically thin discs. Their small line formation radii 
give rise to enhanced response at small time-delays when 
compared to those of H/3 and Mg II. 



APPENDIX B: ANISOTROPIC ILLUMINATION 

A strong assumption of our model is that the gas on the 
bowl-surface is illuminated by a continuum whose shape 
is independent of scale height. Yet, in the TOR model of 
Kawaguchi and Mori (2010,2011) the bowl-shaped geometry 
is formed because the disc emission is strongly anisotropic, 
thereby allowing dust grains to form at much smaller radii 
at low elevations (near to mid-plane of the disc). If an 
anisotropic continuum source is indeed responsible for shap- 
ing the bowl geometry, then it will almost certainly have a 
strong effect on the shape of the 1-d and 2-d responsivity- 
weighted response function and variable emission-line pro- 
file. Here we look at 3 alternative models for the disc con- 
tinuum emission : strong, intermediate and weak anisotropy. 
For each we model the broad band continuum as the sum of 
a varying UV component and a constant X-ray component, 
the relative fractions of which are determined by the depen- 
dence of the intensity of the UV component on the polar 
angle. 

To represent the strong anisotropy dependence we 
adopt the disc-illumination function first proposed by Net- 
zer (1987), ie. 
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Figure Al. The responsivity- weighted 2-d response functions for 
our fiducial model for viewing angles in the range 2-40 degrees. 
TDS and GR are both included, and the turbulence parameter is 
set to frtur-b = 2 - 
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Figure A2. The corresponding variable emission-line profiles 
(solid black line) for ngure lAll Shown in red is the equivalent vari- 
able line profile for emissivity- weighting only (77(7*) = constant) . 
All profiles have been normalised to their peak intensity to aid 
comparison. 



I(4>) = If- cos cf) (1 + 2 cos 0) + J x , (Bl) 

where I uv and I x are the relative intensities (or ionising 
photon fluxes) of the UV and X-ray components, and is 
the polar angle. The weak anisotropy function is taken from 
Nemmen and Brotherton (2010), for which we assume 



for polar angles 



72 degrees, and 



I{4>) = iuv + ix , 



(B2) 



(B3) 



otherwise. Finally, we model an intermediate anisotropic il- 
lumination function with a simple cosine dependence, ie. 



I(<t>) = iuvCOS(/> + / x 



(B4) 



Note we do not calculate new photoionisation model grids 
for a broad range of continuum shapes. Instead as a first ap- 
proximation, we assume that we can use our original emis- 
sivity grids (line flux as a function of radial distance R 
and hydrogen ionising photon flux $h) as a look-up table. 



That is we use the line flux at a given radial distance to 
find the corresponding hydrogen ionising photon flux 
<J>h is then modified according to equations B1-B4, and a 
new flux determined from the same model grid. In figure iBll 
we illustrate the modified H/3 radial emissivity distribution 
for each of the disc illumination functions. Figure IB2I illus- 
trates the 1-d responsivity-weighted response functions (up- 
per panel) and variable emission-line profiles (lower panel) 
resulting from the modified radial emissivity distributions 
as described above for our fiducial model and i = 30 de- 
grees. For each model we assume that the UV continuum 
comprises 99% of the total ionising continuum intensity, ie. 
iuv = 0.99Jtoti. 



APPENDIX C: THE DRIVING CONTINUUM 
LIGHT-CURVE 

The X-ray light curves of AGN show correlated variability 
over a broad range in timescales. This variability is normally 
quantified in terms of the variability power, P, as a function 
of temporal frequency zy, the power spectral density (PSD) 
distribution. For AGN, the X-ray PSD is approximated by 
a power- law in frequency (P oc with slope a = 1 at 
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Figure A3. The corresponding 1-d responsivity- weighted re- 
sponse functions (solid black line) for figure IA1I Shown in red 
is the equivalent emissivity-weighted (ie. 77(7*) = constant) 1-d 
response function. Response functions have been normalised to 
their peak intensities to aid comparison. Responsivity-weighting 
significantly modifies both the form and amplitude of the 1-d re- 
sponse. 



low frequencies, breaking to a slope of 2 at the highest fre- 
quencies (McHardy et al. 2004) appropriate for a red-noise 
process. 

At optical wavelengths, poor temporal sampling has of- 
ten negated the use of the PSD in determining the optical 
continuum variability. Instead, variability is generally char- 
acterised in terms of the first order structure function S(r) 
(Collier and Peterson 2001), where 

1 



S(r) 



(Cl) 



where f(U) is the flux measured at time U, N(r) is the 
number of pairs of points, and the sum is over all pairs for 
which t — tj — 1%. Paltani (1999) showed that the first order 
structure function is related to the one-sided power density 
spectrum P(f) via: 



S(t) 



poo poo 

/ P(f)df- / cos(27r/r)d/ 
Jo Jo 



(C2) 



The form of the structure function can also be approxi- 
mated by a power-law on intermediate timescales, breaking 
to a natter slope on both short timescales (r m in), where 



Figure Bl. The radial surface line emissivity distributions for 
our disc illumination functions: Isotropic continuum - black solid 
line, Netzer's disc illumination function - red dotted line, a disc 
illumination pattern approximating that of Nemmen and Broth- 
erton (2010) - blue dashed line, and a disc illumination function 
with simple cosine dependence - green long-dashed line. 



its value approaches twice the noise variance (cr^), and on 
long-timescales (r ma£E ), where its value approaches twice the 
signal variance (cr^ ar ). 

Collier and Peterson (2001) performed a structure func- 
tion analysis of the UV and optical light-curves of a small 
sample (13) of AGN, including NGC 5548, observed as part 
of the 'AGN Watch' and OSU AGN monitoring programmes. 
They showed that the UV and optical SFs are similar on 
timescales of 5-60 days. In particular, for NGC 5548, the 
UV and optical power-law slope b of the SF on intermediate 
timescales is approximately 1.5, and flattens on a character- 
istic timescale r c har = 40 days. For stationary time series, 
the power-law slope of the structure function b is related 
to the power- law slope of the PSD a, by a = b + 1. Thus, 
for NGC 5548, the slope of the UV and optical SF suggests 
a ~ 2.5. However, we caution that non- stationary effects can 
steepen the derived power-law slope while deviations from a 
strictly power-law slope may be introduced by bias result- 
ing from binning and irregular sampling of the data (Collier 
and Peterson 2001). For a stationary random process S(r) 
is simply related to the autocorrelation function ACF(r), 
such that 



5(r) = 2 [a 2 - ACF{r) 



(C3) 
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Figure B2. 1-d responsivity-weighted response functions and 
variable emission-line profiles for the disc illumination functions 
described in the text for our fiducial bowl-shape BLR geome- 
try, turbulence parameter b tU rb — 2, and observed at inclination 
i = 30 degrees. Isotropic continuum - solid line, Netzer's disc 
illumination function - dotted line, a disc illumination pattern 
approximating that of Nemmen and Brotherton (2010) - dashed 
line (nearly coincident with the solid line), and a disc illumination 
function with a simple cosine dependence - long-dashed line. 



Figure CI. Upper panel - a typical structure function for one of 
our simulated light-curves (open squares). Also shown is the struc- 
ture function derived for the GEX extracted 1337Acontinuum 
from the 1989 IUE monitoring campaign of NGC 5548. Middle- 
panel - the corresponding simulated light-curve sampled at 4 day 
intervals. Lower-panel - the 1337Acontinuum light-curve from the 
1989 monitoring campaign of NGC 5548. 



Thus the form of the SF can vary from one observing sea- 
son to the next even when the process responsible for the 
variability remains the same. Similar critiques have been 
made concerning the use of structure functions in quantify- 
ing Blazar variability (Emmanoulopoulos et al. 2010), and 
in particular when assigning characteristic timescales to fea- 
tures observed in the SF on long timescales where the SF is 
poorly defined. 

In figure IC1I (top panel) we indicate a typical structure 
function calculated from one of our model light-curves (open 
squares). Also shown is the structure function derived from 
the GEX extracted 1337A continuum light-curve (filled tri- 
angles) of NGC 5548 taken from the 1989 IUE monitoring 
campaign (Clavel et al. 1991). The dashed diagonal line in- 
dicates the slope of the structure function over the region of 
interest. In the middle panel we show the simulated light- 
curve sampled at 4 day intervals (to match the sampling 
of the light-curve in the IUE campaign). The lower panel, 
shows the continuum light-curve from the 1989 IUE moni- 
toring campaign. 



